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ABSTRACT 


An  experimental  investigation  of  a  wing  boundary  layer  subjected  to  periodic  tur¬ 
bulent  flow  at  a  Reynolds  number  of  500,000  was  conducted.  Xon-thrusting  turbulence 
pulses  were  generated  at  a  rate  of  50  pulses  per  second  with  a  turbulence  intensity  near 
10%  by  spinning  a  rod  in  the  free  stream  flow  upwind  of  the  test  airfoil. 

Time-varying  velocity  measurements  were  made  at  three  representative  chord  lo¬ 
cations  (laminar,  transitional  turbulent,  fully  turbulent)  at  one  angle  of  attack  using 
single-element  hot-wire  anemometry.  The  characteristics  of  the  boundary  layer  velocity 
profiles,  turbulence  intensity  profiles,  and  velocity  spectra  and  total  spectral  power  were 
documented. 

The  time-varying  boundary  layer  response  could  be  characterized  by  three  flow  re¬ 
gimes:  undisturbed  flow;  turbulence  pulse;  and  recovery  period.  The  boundary  layer 
exhibited  a  cyclic  transition  response  varying  from  the  undisturbed  flow  regime  to  the 
turbulence  pulse  regime  back  to  the  undisturbed  flow  regime.  The  turbulence  pulse  was 
found  to  penetrate  into  the  entire  boundary  layer.  The  turbulence  pulse  prompted 
laminarization  of  the  flow  during  the  recovery  period.  Laminarization  apparently  re¬ 
sulted  from  rapid  acceleration  of  the  near-surface  flow  within  the  boundary  layer  due  to 
turbulence  induced  momentum  transfer  down  through  the  boundary  layer  and  local  flow 
acceleration  following  the  velocity  deficit  of  the  turbulence  pulse.  The  effects  of  the 
periodic  turbulence  pulse  were  most  noticeable  in  the  recovery  period  at  the 
transitional  turbulent  and  fully  turbulent  regions  of  the  wing.  Recovery  was  generally 
characterized  by  lowered  turbulence  levels  and  reduction  in  boundary  layer  thickness, 
indicating  a  stabilizing  effect  of  the  external  turbulence  pulse. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Airfoil  experimental  research  at  low  chord  Reynolds  numbers  has  become  a  major 
field  of  academic,  military,  and  civil  interest  in  the  last  few  years  due  to  the  current  de¬ 
velopment  of  low-speed  remotely  piloted  vehicles  (RPVs),  manned  and  unmanned 
ultra-light  aircraft,  human-powered  aircraft,  and  high  altitude  station  keeping  aircraft. 
Low  Reynolds  number  flight  (Rc  <  500,000)  is  generally  considered  as  very  low  speed 
flight,  or  flight  with  a  short  wing  chord.  Low  Reynolds  number  flight  also  occurs  at  high 
altitudes  where  the  air  density  is  low.  Most  low  Reynolds  number  data  in  the  past  have 
come  from  experiments  on  gliders  and  model  airplanes.  [Ref.  1] 

The  aerodynamic  response  of  a  wing  at  low  Reynolds  numbers  is  considerably  dif¬ 
ferent  than  for  Reynolds  numbers  greater  than  about  500.000.  At  low  Reynolds  num¬ 
bers.  laminar  flow  exists  for  a  distance  downstream  from  the  leading  edge.  As  such,  the 
laminar  to  turbulent  transition  occurs  much  further  back  on  the  wing,  often  in  the  ad¬ 
verse  pressure  gradient  region.  The  transition  area  is  often  quite  broad  and  not  well 
defined.  In  addition,  the  transition  process  is  very  sensitive  to  free  stream  disturbances, 
external  turbulence,  and  surface  imperfections. 

Laminar  (low  tends  to  separate  easily  at  low  Reynolds  numbers  especially  in  the 
presence  of  an  adverse  pressure  gradient.  Leading  and  trailing  edge  separation  is  fre¬ 
quently  observed  as  well  as  separation  along  the  airfoil.  Separation  along  the  airfoil  is 
commonly  referred  to  as  the  "laminar  separation  bubble".  Laminar  separation  bubbles 
greatly  affect  the  performance  of  the  wing.  The  most  efficient  low  Reynolds  number 
airfoils  are  designed  to  obtain  maximum  benefit  from  the  separation  bubble;  but  only 
because  the  separation  bubble  cannot  be  completely  eliminated  at  low  Reynolds  num¬ 
bers.  Numerous  factors  also  affect  the  size,  shape,  and  character  of  the  separation 
bubble  including  free  stream  disturbances,  vibration,  surface  condition,  angle  of  attack, 
etc. 

The  critical  nature  of  the  boundary  layer  and  the  many  factors  that  influence  it  have 
prompted  many  recent  investigations  in  the  characteristics  of  airfoils  and  wings  at  low 
Reynolds  numbers.  Although  investigations  under  steady  state  conditions  provide  im¬ 
portant  lift,  drag,  and  airfoil  performance  data,  they  may  not  reliably  predict  the  laminar 
to  turbulent  transition,  leading  and  trailing  edge  separation,  and  formation  and  reat- 
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tachment  of  the  separation  bubble  as  these  phenomena  may  be  influenced  to  a  high  de¬ 
gree  by  unsteadiness  in  the  (low.  [Ref.  1] 

Unsteady  low  Reynolds  number  aerodynamic  research  has  recently  been  undertaken 
on  three  fronts: 

Airfoil  response  to  high  steady-state  free  stream  turbulence 

Airfoil  response  to  low  to  moderate  frequency  sinusoidal  modulation  of  the  free 
stream 

Airfoil  response  to  transient  free  stream  disturbances 

Meier  and  Kreplin  [Ref.  2],  Potter  et  al.  [Ref.  3],  and  Kindelspire  [Ref.  4].  among 
others,  have  investigated  airfoil  performance  and  boundary  layer  development  and  be¬ 
havior  in  the  presence  of  high  free  stream  turbulence.  Meier  and  Kreplin  found  that 
both  the  intensity  and  spectral  distribution  (length  scale)  of  the  free  stream  turbulence 
affected  airfoil  performance.  The  most  significant  changes  were  noted  when  the  length 
scales  were  on  the  order  of  the  boundary  layer  thickness.  Potter  et  al.  also  found  that 
the  free  stream  turbulence  intensity  spread  into  the  entire  boundary  layer,  but  the  effect 
was  dependent  on  length  scale.  They  also  noted  that  the  free  stream  turbulence  has  to 
penetrate  and  alter  the  velocity  profile  before  producing  any  significant  change  in  the 
boundary  layer  detachment  location.  Kindelspire  found  that  increasing  levels  of  free 
stream  turbulence  caused  correspondingly  earlier  laminar  to  turbulent  transition,  lie 
also  confirmed  that  urbulence  length  scales  on  the  order  of  the  boundary  layer  thickness 
had  the  most  influence  on  boundary  layer  thickness. 

Brendel  [Ref.  5]  performed  a  wind  tunnel  study  with  unsteady  flow  consisting  of  a 
streamwise  sinusoidal  velocity  variation.  Brendel  found  that,  in  general,  the  imposed 
unsteadiness  had  little  effect  on  the  mean  boundary  layer  thickness.  The  unsteadiness 
did,  however,  reduce  the  thickness  and  length  of  the  separation  region  associated  with 
the  separation  bubble.  He  also  found  that  lower  frequencies  (large  length  scales)  had 
no  effect  on  the  overall  flow. 

Howard  [Ref.  6[  most  recently  investigated  the  time  vary  ing  characteristics  of  the 
boundary  layer  of  an  airfoil  immersed  in  a  propeller  slipstream  (see  Figure  1).  The 
propeller  wake  subjected  the  airfoil  to  a  transient  free  stream  disturbance.  Howard 
found  that  the  wing  boundary  layer  varied  from  laminar  to  turbulent  in  response  to  the 
passage  of  the  propeller  wake.  He  also  found,  in  general,  that  the  effects  of  the  propeller 
slipstream  tended  to  enhance  the  stability  of  the  boundary  layer  and  reduced  the  drag 
coefficient  in  the  laminar  portion  of  the  slipstream  cycle  below  its  undisturbed  value. 
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Laminar 


Figure  1.  Wing  Immersed  in  Propeller  Slipstream:  l Ref.  6:  p.  17) 


In  particular,  Howard  noted  a  distinct  relaminarization  process  at  work  during  the 
transition  from  turbulent  to  laminar  flow  as  evidenced  by  a  decrease  in  turbulence  and 
an  increase  in  velocity  near  the  surface. 

B.  THESIS  OBJECTIVES 

The  goal  of  this  study  is  to  thoroughly  investigate  the  boundary  layer  response  over 
a  wing  to  periodic  free  stream  disturbance  at  the  laminar,  transitional,  and  turbulent 
portions  of  the  wing. 

This  thesis  describes  the  design,  fabrication,  and  testing  of  a  non-thrusting  periodic 
turbulence  generating  device  tc  simulate  transient  phenomena  such  as  wind  gust  or  tur¬ 
bulent  wake  disturbances.  The  investigation  consists  of  wind  tunnel  tests  of  an  airfoil 
operating  in  a  periodically  disturbed  free  stream  at  a  chord  Reynolds  number  of  ap¬ 
proximately  500,000.  Measured  data  are  presented  in  the  form  of  non-stationary  mean 
velocity  and  turbulence  time  histories  at  various  heights  above  the  surface;  as 
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ensemble-average  velocity  profiles  and  turbulence  profiles:  and  as  velocity  spectra  and 
total  power  plots. 

These  data  will  form  a  base  for  follow-on  research  at  high  angles  of  attack  and  for 
work  within  the  separation  region  of  low  Reynolds  number  airfoils. 


II.  EXPERIMENTAL  EQUIPMENT 


A.  APPARATUS 

The  major  hardware  items  utilized  in  this  investigation  included  the  wind  tunnel,  the 
wing  section  and  wing  mount  assembly,  the  periodic  turbulence  generator,  the  hot-wire 
anemometer,  the  traversing  mechanism,  and  the  data  acquisition  and  processing  system. 
This  chapter  presents  the  specifications,  design,  setup,  and  operating  characteristics  of 
all  equipment  used. 

1.  Wind  Tunnel 

The  Naval  Postgraduate  School  low  speed  horizontal  wind  tunnel  was  employed 
for  this  experiment.  Figure  2  shows  the  tunnel  layout. 

The  tunnel  is  of  single  return  design  and  measures  64  feet  in  length,  25.5  feet  in 
width,  and  supports  a  32-inch  by  45-inch  by  4-foot  test  section.  The  tunnel  was  designed 
by  Aerolab  Development  Company  of  Pasadena,  California  and  has  been  in  continuous 
use  at  the  Navy  School  since  the  middle  1950's.  [Ref.  7] 

The  tunnel  consists  of  five  basic  units:  the  fan  section:  diffusers;  settling  cham¬ 
ber;  contraction  cone;  and  test  section.  The  tunnel  is  constructed  of  wood  throughout 
and  recently  received  a  coat  of  fire  resistant  paint  on  the  interior  surface  to  meet  OSH  A 
standards.  The  paint  roughened  the  surface  somewhat:  however,  the  level  of  free-stream 
turbulence  was  not  significantly  worsened  due  to  the  numerous  anti-turbulence  features 
built  into  the  tunnel. 

The  tunnel  is  powered  by  a  100-horsepower  electric  motor  driving  a  three-bladed 
variable-pitch  propeller.  Tunnel  velocity  is  controlled  by  electrically  setting  blade  pitch 
angle  at  the  control  unit  located  adjacent  to  the  test  section.  A  four-speed  Dodge 
transmission  is  installed  between  the  motor  and  the  fan  to  provide  greater  tunnel  velocity 
control  over  the  design  range  of  the  tunnel  (up  to  290  ft  sec).  The  transmission  was  left 
in  3rd  gear  for  all  data  runs  with  a  fan  shaft  speed  of  750  RPM. 

Eight  radial  stator  blades  are  mounted  directly  downstream  of  the  fan  on  the 
drive  shaft  housing.  These  stator  blades  serve  to  straighten  the  fiow  and  remove  the 
swirl  generated  by  the  fan.  These  arc  the  first  of  the  installed  anti-turbulence  devices. 

The  diffuser  sections  are  located  on  either  side  of  the  fan  section.  The  gradual 
increase  in  cross-sectional  area  and  the  long  length  of  the  diffuser  sections  gradually  slow 
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Section  B-B  Section  A-A 

Figure  2.  NPS  Lou  Speed  Wind  Tunnel 


the  air  while  introducing  little  additional  turbulence  and  keeping  friction  losses  to  a 
minimum. 

The  large  settling  chamber,  located  at  the  left  of  Figure  2,  damps  any  remaining 
flow  unsteadiness.  In  addition,  two  fine  wire  screens  are  mounted  six  inches  apart  in  the 
settling  chamber.  These  serve  to  further  break  up  turbulent  fluctuations  in  the  flow 
prior  to  the  air  entering  the  test  section. 

The  contraction  cone  accelerates  the  air  and  delivers  it  to  the  test  section.  The 
large  contraction  ration  of  10:1  helps  assure  low  turbulence  in  the  test  section. 

The  test  section  of  the  low  speed  wind  tunnel  is  equipped  with  a  full  size  re¬ 
flection  plane  mounted  4  inches  above  the  tunnel  floor.  The  reflection  plane  is  perma¬ 
nently  mounted  and  has  an  electrically-driven  turntable  mechanism  in  the  center  (not 
used  in  this  investigation).  The  test  section  is  additionally  equipped  with  full  length 
lighting  located  in  the  corner  fillets  and  has  a  hinged  window  on  either  side.  The  walls 
of  the  test  section  are  slightly  divergent  to  counteract  the  contraction  resulting  from 
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boundary  layer  growth.  Breather  slots  are  installed  immediately  downstream  of  the  test 
section  to  maintain  the  test  section  at  atmospheric  pressure  and  to  allow  the  tunnel  to 
breath  as  necessary  on  start  up  and  shutdown. 

Two  sets  of  static  pressure  port  rings  are  installed  in  the  tunnel  and  used  to 
measure  the  tunnel  speed.  One  ring  set  is  installed  just  upstream  of  the  contraction  cone 
and  the  other  set  is  installed  at  the  test  section  entrance.  The  pressure  differential  across 
the  two  ring  sets  is  measured  on  a  water  micro-manometer.  This  pressure  differential, 
once  calibrated,  corresponds  to  the  velocity  in  the  test  section.  The  wind  tunnel  velocity 
calibration  procedure  is  covered  in  Chapter  III. 

2.  Wing  Section  and  Wing  Mount  Assembly 

The  wing  section  and  wing  mount  assembly  are  shown  in  Figure  3.  This  in¬ 
stallation  was  designed  by  Kindelspire  [Ref.  4:  pp.  11-15]  and  a  complete  description  of 
the  apparatus  is  contained  in  his  work.  The  following  paragraphs  highlight  the  pertinent 
aspects  of  the  wing  and  wing  mount  system. 

The  airfoil  tested  was  a  section  of  a  Bell  helicopter  tail  rotor  (model  and  type 
unknown)  having  a  10-inch  chord  and  24-inch  span.  The  airfoil  shape  is  shown  in  Fig¬ 
ure  4  and  the  section  coordinates  are  given  in  Appendix  A. 

The  airfoil  section  was  stripped  to  bare  metal,  filled,  sanded,  and  finished  with 
high  gloss  black  lacquer  prior  to  installation.  The  lacquer  finish  was  polished  using  fine 
polishing  compound  and  given  a  final  coat  of  automobile  wax  to  bring  it  to  the  highest 
possible  gloss.  A  glossy,  reflective  surface  was  important  as  the  hot-wire  anemometer 
probe  was  positioned  to  within  0.0032  inches  of  the  surface  by  viewing  the  probe  and  its 
reflection  through  a  sighting  level.  The  probe  positioning  procedure  is  discussed  in 
Chapter  III. 

The  wing  mount  assembly  consists  of  two  8-in  wide  machined  aluminum 
endplates  mounted  vertically  between  the  reflection  plane  and  the  top  of  the  test  section. 
The  wing  section  is  mounted  between  the  endplates  on  rotatable  inserts  that  allow  the 
angle  of  attack  (AOA)  to  be  varied  from  +35  degrees  to  -15  degrees.  Graduated 
markings  on  the  rotatable  inserts  are  used  for  coarse  setting  of  AOA.  The  acutal  AOA 
is  determined  by  measuring  the  height  of  the  leading  and  trailing  edges  above  the  re¬ 
flection  plane  and  solving  analytically  for  AOA, 

All  data  collected  during  this  experiment  were  taken  l-inch  left  of  mid-span 

% 

(looking  from  the  rear)  due  to  the  presence  of  several  mars  in  the  wing  surface  at  the 
mid-span  line  which  were  caused  by  the  measuring  micrometer  when  the  wing  section 
coordinates  were  documented.  The  1-inch  offset  was  sufficient  to  remain  free  of  dis- 
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Figure  4.  Experimental  Airfoil  Section 

The  support  assembly  was  constructed  out  of  aluminum  stock  and  designed  to 
bolt  to  the  reflection  plane  such  that  the  hub  was  position  2.5  inches  below  and  1  inch 
ahead  of  the  wing  at  center  span.  The  hub  was  placed  well  below  the  wing  in  order  to 
allow  the  disturbances  generated  by  the  hub  section  to  pass  under  the  wing.  The  leading 
edges  of  the  mount  were  also  rounded  and,  or  beveled  to  minimize  drag  and  free  stream 
disturbances.  Figure  6  shows  the  wing  wing  mount  and  periodic  turbulence  generator 
installed  in  the  NPS  wind  tunnel. 

The  vertical  strut  of  the  support  assembly  was  bored  to  accept  the  output  end 
of  the  flexible  drive  shaft.  The  drive  shaft  was  held  firmly  in  place  by  a  bracket  mounted 
on  the  backside  of  the  vertical  strut. 

The  bar  and  hub  were  both  machined  out  of  stainless  steel.  The  18-inch  by 
3  8-inch  diameter  bar  was  secured  to  the  hub  by  a  standard  1/8-inch  roll  pin.  The  hub 
was  bored  and  tapped  to  screw  onto  the  threaded  end  of  the  flexible  drive  shaft  extend¬ 
ing  from  the  support  assembly.  The  hub  was  secured  in  place  on  the  threaded  shaft  by 
a  locknut.  In  addition,  the  counterclockwise  direction  of  rotation  was  chosen  such  that 
the  hub  would  tend  to  tighten  during  operation. 

The  38-inch  flexible  drive  shaft  used  to  couple  the  motor  and  the  bar/hub  was 
purchased  at  a  local  hardware  store  as  an  electric  drill  extension.  As  noted  previously, 
the  output  end  had  a  3/S-inch  threaded  shaft  onto  which  the  hub  assembly  was  mounted. 
The  input  end  of  the  flexible  shaft  was  made  to  fit  a  1 /4-inch  drill  chuck.  Use  of  the 
flexible  drive  shaft,  rather  than  direct  drive,  reduced  the  size,  weight,  and  complexity  of 
the  support  assembly  as  it  did  not  have  to  support  the  heavy  drive  motor. 

The  drive  motor  and  position  sensor  are  shown  in  Figure  7.  The  motor  was 
mounted  on  the  tunnel  floor  well  clear  of  the  wing  environment  at  the  downwind  end 
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Periodic  Turbulence  Generating  Assembly 


Figure 


ion.  I  he  motor  is  a  1  4-horsepower  permanent  magnet  variable  speed 
'need  could  he  controlled  over  the  ranee  wt-25'io  RPM.  An  adapter  with 


.vas  manufactured  to  lit  the  1 


ft  speed  and  shaft  position  from  a  toothed  steel  d>K  mounted  01 
octhcd  steel  disk,  seen  in  1  mure  was  modified  bv  removing 


of  tite  original  teeth.  1  iuis  modified,  an  electrical  pulse  was  generated  tw  icc  c  ti  c  n 
revolution  corresponding  to  each  passage  of  the  bar  in  front  of  the  wing  section.  ‘I  his 
signal  was  used  to  drive  an  electronic  RPM  counter  and  pros  idea  timing  information  for 
the  data  acquisition  system  and  data  analysis  software. 

4.  Hot-Wire  Anemometer 


Velocitv  and  turbulence  measurements  were  made  using  a  I  hermo  System.  Inc 
t  f  S I  >  II  A-l""  S\  stem  Mode!  1 5* »  (  onstunt  Temperature  Anemometer  and  Model  f 


Figure  6.  Wing/Wing  Mount  and  Periodic  Turbulence  Generator  Installed  in  NFS 
Wind  Tunnel 


Sign;;!  Conditioner.  Figure  S  shows  the  I  FA- 1  no  system  installation.  The  ITA-1"" 
s\  stern  is  used  with  touch  pane!  entry  el' all  anemometer  and  signal  conditioner  settings 
ane  digit..!  panel  readouts. 

The  Model  IS"  Ar.emoirieter  module  was  operated  in  the  standard  bridge  con¬ 
figuration  with  a  standard  15-foot  RG-5SA  l  probe  cable.  This  configuration  prenided 
extremely  low  noise  output  signals  with  a  typical  frequency  response  of  DC  to  15"  KI !/ 
[Ref.  S:  pp.  S-?  to  S — » j.  These  parameters  were  important  in  order  to  accurately  measure 
the  low  level,  high  frequency  turbulence  components  of  the  velocity  signal. 

Tire  hot-wire  signal  was  processed  using  the  Model  157  Signal  Conditioner  to 
improve  the  veioeitx  resolution  of  the  measurement  system  and  to  provide  low  pass  fil¬ 
terin'.:.  live  unconditioned  hot-wire  output  voltage  oxer  the  calibration  \ elocity  range 
varied  from  1.5  volts  at  a  flow  rate  of  4"  ft  sec  to  l.S  xolts  at  14"  ft  see.  With  the  g:\en 
anal.  _  to  cheital  xcdtaee  resolution  of  2.44  millivolts,  the  resulting  \ elocity  resolution 
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Figure  7.  Drive  Motor  and  Position  Sensor 

was.  at  best.  O.S  ft  sec.  By  applying  a  -1.0  volt  offset  and  gain  of  5  to  the  raw  signal  with 
the  signal  conditioner,  the  output  voltage  over  the  calibration  range  was  made  to  vary 
from  2.5  to  4.0  volts  giving  a  velocity  resolution  of  0.16  ft  sec-thc  improvement  in  re¬ 
solution  being  directly  related  to  the  gain  setting. 

The  low  pass  filter  of  the  Model  157  Signal  Conditioner  was  set  as  close  as 
possible  to  the  Nyquist  cutoff  frequency.  The  Nyquist  cutoff  frequency.  fn  ,  is  defined 
by  the  formula: 


where  f,  is  the  digital  sampling  rate  in  Hertz  (Hz).  To  avoid  serious  errors  in  the  data, 
it  is  important  that  the  signal  have  no  significant  components  in  the  frequency  range 
above  fn  prior  to  digital  conversion  [Ref.  9:  pp.  68-69].  Setting  the  low  pass  filter  to  fn 
effectively  removes  all  the  unwanted  high  frequency  components. 
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Figure  8.  IFA- 100  Anemometer  System 

A  Dantec  55p  1 1  straight  general  purpose  hot-wire  probe  with  a  5-micron  di¬ 
ameter  platinum  plated  tungsten  wire  was  used  for  all  measurements. 

The  oscilliscope  and  digital  voltmeter,  also  mounted  on  the  1FA-100  equipment 
rack  (Figure  8),  were  used  to  monitor  the  hot-wire  signal  during  each  run.  The  digital 
voltmeter  was  set  to  read  the  mean  DC,  or  steady  state,  component  of  the  signal  and 
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served  to  monitor  the  steady  state  velocity.  The  oscilliscope  displayed  the  AC.  or  tur¬ 
bulent,  component  of  the  signal  and  served  to  monitor  turbulence  in  the  How. 

The  calibration  procedures  for  the  anemometer  are  covered  in  Chapter  Ill  of 
this  report. 

5.  Three  Dimensional  Traverser 

Figure  9  shows  the  Velmex  8300  Control  Driver  3-D  traversing  system  em¬ 
ployed  to  position  the  hot-wire  probe  through  the  boundary  layer.  The  traverser  con¬ 
figuration  used  by  this  investigator  was  identical  to  that  utilized  and  described  by 
Kindelspire  [Ref.  4:  pp.  21-24]. 

The  traverser  was  mounted  above  the  test  section,  thus  allowing  the  hot-wire 
probe  to  be  positioned  spanwise.  chordwise,  and  vertically  in  relation  to  the  wing  sec¬ 
tion.  Minimum  incremental  movement  in  any  of  the  three  directions  was  0.000125 
inches.  The  traverser  could  be  controlled  either  manually  or  by  software  command  from 
a  personal  computer  interface.  An  IBM  PC  AT  was  used  to  control  the  traverser  in  this 
installation.  The  software  routine  used  to  drive  the  traverser  is  explained  in  a  later  sec¬ 
tion  of  this  chapter. 

6.  Analog-to-Digital  Data  System 

Data  collection  was  accomplished  using  a  PC  AT-type  personal  computer  cou¬ 
pled  to  a  Metrabyte  Model  DAS-16F  analog-to-digital  converter. 

The  PC  AT-type  personal  computer  used  for  data  collection  and  analysis  was 
configured  with  a  40  megabyte  hard  disk,  two  high  density  disk  drives,  and  operated  at 
a  throughput  speed  of  12  MHz.  The  computer  was  additionally  equipped  with  1  mega¬ 
byte  of  random  access  memory.  The  large  storage  capacity  was  necessary  to  handle  the 
volume  of  data  collected  in  this  investigation.  Each  sample  point  yielded  close  to  1 
megabyte  of  raw  data  and  over  150  points  were  sampled. 

The  Metrabyte  DAS-16F  is  a  12-bit  very  high  speed  analog  to  digital  interface 
designed  for  use  with  IBM  and  IBM  compatible  personal  computers.  The  board  fits 
directly  in  an  expansion  slot  in  the  computer.  The  analog  inputs  were  brought  in 
through  an  external  interface  junction  box.  The  board  operates  at  a  maximum  rate  of 
100,000  samples  sec  (Hz)  with  up  to  16  single  ended  or  8  differential  analog  input 
channels.  The  DAS-16F  offers  a  wide  range  of  operating  parameters  including  four 
triggering  options,  three  data  transfer  options,  nine  selectable  input  voltage  ranges,  and 
various  other  input  output  modes. 

The  DAS-16F  was  configured  to  read  two  differential  channels  (hot-wire  voltage 
and  timing  pulse)  over  the  input  range  of  0-5  volts.  This  resulted  in  an  analog-to-digital 
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Figure  9.  3-D  Traverser 

resolution  of  2.44  millivolts.  The  board  was  operated  at  its  maximum  sample  rate  of 
1 1 1(1.000  \  \z  (So.oon  11/  per  channel)  during  this  experiment. 

B.  SOFTWARE 

The  computer  played  an  important  role  in  all  phases  of  data  collection  and  analysis 
in  this  investigation.  .  The  entire  experiment  was  conducted  using  two  PC  AT-type 
computers.  The  software  programs  created  for  each  phase  of  the  investigation  are  de¬ 
tailed  in  this  section.  The  computer  code  for  each  of  the  programs  developed  by  this 
author  is  contained  in  Appendix  C  through  Appendix  F. 

1.  Data  Collection  Software 

l  he  programs  in  this  section  facilitated  data  collection  and  performed  the  fol¬ 
lowing  control  functions: 

3-D  Traverser  Control  (TRAVERSE. BAS) 

I  Iot-Wire  Calibration  (HWCAE.BAS) 
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Analog-to-Digital  Board  control  (LABTECH  NOTEBOOK) 

Data  Acquisition  in  the  Boundary  Layer  (BLDATA. BAS) 

a.  Program  TRAVERSE. BAS 

The  program  TRAVERSE. BAS  provides  discrete  or  automated  motor 
movement  control  of  the  YELMEX  8300  traverser  through  interactive  keyboard  com¬ 
mands  from  an  IBM  PC  AT  computer.  The  TRAVERSE. BAS  program  written  by 
Kindelspire  [Ref.  4:  pp.  27-30J  was  modified  by  this  author  to  remove  manual  data  entry 
routines  and  to  improve  the  accuracy  of  the  motor  stepping  routines. 

Program  TRAVERSE  gives  the  operator  the  option  to  input  individual 
discrete  step  commands  from  the  keyboard  or  to  elect  automated  movement  sequences. 
Discrete  inputs  were  generally  used  to  position  the  hot-wire  in  the  boundary  layer  prior 
to  each  investigation.  Discrete  movements  could  be  entered  for  any  or  all  of  the  three 
directions  with  the  minimum  movement  being  0.000125  inches  in  any  one  direction. 
Once  the  hot-wire  was  positioned,  the  operator  then  selected  computer  controlled 
movement  to  advance  the  probe  through  the  boundary  layer  during  the  data  acquisition 
process. 

The  computer  stepped  the  hot-wire  through  the  boundary  layer  according 
to  the  positioning  algorithm  of  Equation  2: 


z(j)  =  BL 


[sin(f )] 


1  + 


(i-D  11 
(n-l)  Jj 


7.(j )  =  Position  in  the  boundary  layer 
BL  =  Estimated  boundary  layer  thickness 
n  =  Number  of  data  points  desired 
j  =  Incremental  data  point 


(2) 


The  sine  function  based  algorithm  of  Equation  2  is  designed  to  concentrate 
many  smaller  incremental  steps  at  the  bottom  of  the  boundary  layer  and  gradually  in¬ 
crease  the  increment  size  toward  the  top  of  the  layer.  Special  care  was  taken  in  writing 
the  actual  stepping  routine  to  ensure  motor  movement  commands  were  given  in  multi¬ 
ples  of  the  motor  step  interval  (0.000125  inches)  to  preclude  cumulative  roundoff  errors 
throughout  the  boundary  layer. 

b.  Program  HWCAL.BAS 

This  program,  written  by  the  author,  is  used  to  calibrate  the  hot-wire 
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anemometer.  The  program  was  written  in  Microsoft 2  QuickBasic  and  compiled  to  gen¬ 
erate  a  stand  alone  executable  program. 

The  hot-wire  is  calibrated  against  a  reference  pitot-static  probe  located  in 
the  test  section  free  stream  adjacent  to  the  hot-wire.  Pitot-static  input  is  in  centimeters 
(cm)  of  water  as  read  on  a  micro-manometer.  These  data  are  manually  entered  on  the 
keyboard  by  the  operator  for  each  calibration  point  while  the  corresponding  hot-wire 
voltage  from  the  I  FA- 100  Anemometer  is  digitally  read  by  the  Metrabyte  DAS-16F 
analog-to-digital  converter  board.  The  operator  also  enters  the  ambient  temperature 
and  pressure  at  the  conclusion  of  each  run.  One  thousand  hot-wire  samples  are  read  and 
averaged  at  each  calibration  point3.  The  hot-wire  was  generally  calibrated  at  1,2,  3,  4, 
6.  8.  and  10  cm  of  water  providing  accurate  calibration  data  over  the  velocity  range 
40-140  ft  sec.  Program  output  consists  of  five  columns  of  data  suitable  for  plotting  us¬ 
ing  most  current  graphing  software  packages: 

Pitot-static  pressure  (cm  of  water) 

Hot-wire  voltage  (volts) 

Frce-stream  velocity  (ft  sec) 

(Hot-wire  voltage)  **  2 
Velocity  **  .5 

The  hot-wire  calibration  procedures  are  covered  in  Chapter  Ill  of  this  report, 
c.  LAB  TECH  NO  T EBOOK 

The  commercial  software  package  LABTECH  NOTEBOOK  f Version  4.3 j 
marketed  by  Laboratory  Technologies  Corporation^  provided  set  up  and  operating 
control  of  the  DAS-16F  analog-to-digital  board.  Using  the  NOTEBOOK  menus,  the 
operator  was  able  to  enter  all  the  operating  parameters  including  sample  rate,  channels 
to  sample,  run  duration,  start  mode,  output  data  file  name  and  output  device,  and  data 
output  format.  NOTEBOOK  allowed  this  investigator  to  take  and  store  up  to  130,000 
data  samples  during  any  one  run  while  operating  at  the  100,000  Hz  rate.  Other  software 
control  programs  appeared  to  be  limited  to  32,000  samples  under  the  same  operating 
conditions.  This  feature  was  paramount,  as  a  minimum  data  capability  of  92.000  dis¬ 
crete  samples  was  required  at  each  point. 

2  Microsoft  is  a  registered  trademark  of  Microsoft  Corporation 

3  The  suggested  analog  to  digital  sample  rate  is  1000  Hz  with  a  run  period  of  1.01  seconds. 

4  LABTFC1I  is  a  registered  trademark  of  Laboratory  Technologies  Corporation. 
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d.  Program  BLDA  TA.BAS 

This  driver  program  controls  the  analog-to-digital  data  acquisition  process. 
Program  BLDATA.BAS  was  written  in  QuickBasic  and  compiled  to  generate  a  stand 
alone  executable  program.  The  user  inputs  the  number  of  data  points  desired  in  the 
boundary  layer  and  a  common  output  file  name.  BLDATA  then  calls  NOTEBOOK 
which  in  turn  initiates  data  collection.  Once  the  data  are  collected  and  the  output  file 
written  by  NOTEBOOK,  BLDATA  renames  the  output  file  giving  it  a  sequentially 
numbered  tag  and  prepares  for  another  data  point  in  the  boundary  layer.  At  the  com¬ 
pletion  of  the  run,  the  investigator  has  a  directory  of  sequentially  named  files  relating  to 
each  data  point. 

2.  Data  Analysis  Software 

The  author  used  the  data  analysis  software  routines  specifically  developed  by 
Johnson  [Ref.  10]  to  analyze  the  periodic,  non-steady  velocity  and  turbulence  data  ob¬ 
tained  in  this  investigation.  These  programs  are  not  reproduced  here  due  to  their  length 
and  considerable  documentation.  A  summary,  however,  of  Johnson's  main  programs 
as  well  as  those  written  by  the  author  is  provided. 

a.  Program  CONVERT. FOR 

This  program  was  written  by  the  author  to  convert  the  raw  hot-wire  voltage 
and  position  sensor  voltage  into  velocity  data  and  a  synchronization  pulse.  Input  data 
was  in  ASCII  integer  form  with  the  integer  values  0  to  4095  representing  0  to  5  volts. 
The  bar  sensor  voltage  was  processed  to  return  a  2  at  the  exact  center  of  the  position 
sensor  pulse  and  a  0  elsewhere.  The  hot-wire  voltage  to  velocity  conversion  constants 
are  obtained  from  the  hot-wire  calibration  data  from  program  HWCAL  discussed  pre¬ 
viously. 

b.  Program  MEA  NS  MU.  FOR 

This  program,  written  by  Johnson  [Ref.  10:  pp.  158-166].  computes  the  time 
varying  mean  of  a  non-stationary,  periodic  time  series.  The  data  are  first  ensemble  av¬ 
eraged  in  the  time  domain  and  then  transformed  into  the  frequency  domain  via  a  Fast 
Fourier  Transform  (FFTt.  The  ensemble  is  then  smoothed  by  applying  a  low  pass  dig¬ 
ital  filter  to  the  ensemble.  The  ensemble  is  then  transformed  back  into  the  time  domain 
by  an  inverse  FFT  and  written  to  the  output  file.  The  amount  of  smoothing  is  a  func¬ 
tion  of  the  low  pass  filter  cutoff  frequency.  The  lower  the  cutoff  frequency,  the  greater 
the  smoothing  and  vice  versa.  The  amount  of  smoothing  can  be  specified  by  the  user 
through  keyboard  input.  The  methodology  of  selecting  the  low  pass  cutoff  frequency  is 
detailed  in  Chapter  IV. 
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c.  Program  T UR  BIX.  FOR 

Written  by  Johnson  [Ref.  10:  pp.  173-181],  this  program  computes  the  in¬ 
stantaneous  turbulence  intensity  of  a  non-stationary,  periodic  time  series.  It  is  the  tur¬ 
bulence  intensity  counterpart  of  the  non-stationary  mean  program  MEANSMU. 
TURBIN  first  calculates  the  time  varying  mean  velocity  of  each  ensemble  by  low  pass 
filtering  in  the  frequency  domain.  TURBIN  then  determines  the  instantaneous  turbulent 
(fluctuating)  component  at  each  point  on  each  ensemble.  The  turbulence  components 
are  then  ensemble  averaged  and  and  written  to  the  output  file.  As  in  MEANSMU,  the 
user  can  specify  the  cutoff  frequency  of  the  low  pass  filter  used  to  smooth  the  individual 
ensembles. 

d.  Program  PSD. FOR 

This  program,  written  by  Johnson  [Ref  10:  pp.  1S2-190],  computes  the 
power  spectral  density  of  a  real  discrete  steady  time  series.  This  program  was  for  spec¬ 
tral  analysis  of  the  free  stream  flow  in  the  test  section. 

e.  Program  EJS'SPSDA  V.FOR 

This  program,  written  by  Johnson  [Ref  10:  pp.  191-202],  performs  spectral 
analysis  on  a  short  time  segment  of  an  ensemble  of  non-stationary.  periodic,  turbulent 
velocity  data.  Whereas  program  PSD  performs  spectral  analysis  on  a  steady  state  signal. 
ENSPSDAY  investigates  the  time  varying  spectrum  on  selected  portions  of  the  ensemble 
itself  As  such,  the  author  was  able  to  separately  investigate  the  undisturbed  flow  seg¬ 
ment.  the  turbulence  pulse,  and  the  laminar  flow  recovery  period  and  characterize  the 
power  spectrum  and  power  content  of  each. 
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III.  EXPERIMENTAL  PROCEDURE 


A.  WIND  TUNNEL  TEST  SECTION  CALIBRATION 

The  purpose  of  wind  tunnel  test  section  calibration  was  to  determine  the  relationship 
between  the  actual  free  stream  dynamic  pressure  (q)  in  the  test  section  and  the  measured 
static  pressure  differential  (Ap)  across  the  contraction  cone.  This  relationship  is  used  to 
indirectly  determine  test  section  free  stream  velocity  by  measuring  Ap. 

For  this  calibration,  a  reference  pitot-static  probe  was  inserted  in  the  center  of  the 
test  section.  A  differential  micro-manometer  was  connected  across  the  total  and  static 
ports  of  the  pitot-static  probe  to  measure  q.  The  relationship  between  the  pitot-static 
probe  differential  pressure,  q,  and  test  section  velocity  is  given  by  Equation  3: 

pv2 

q  =  (Pt  -  p)  =  —  (3) 

p,  =  Total  pressure  measured  by  the  pitot  tube 
p  =  Static  pressure  in  the  test  section 
q  =  Test  section  dynamic  pressure 
p  =  Density  of  air 
v  =  Velocity  in  the  test  section 

The  static  pressure  differential  across  the  contraction  cone,  Ap,  was  measured  by  a 
second  differential  micro-manometer  connected  across  the  two  static  port  ring  assem¬ 
blies  mounted  in  the  wind  tunnel  walls  just  before  and  after  the  contraction  cone.  In 
most  cases  q  is  linearly  related  to  Ap  and  can  be  written: 

q  =  F(Ap)  (4) 

where  F  is  a  constant  called  the  Tunnel  Calibration  Factor  [Ref.  7:  pp.  Al-1  to  A 1-5]. 
F  is  determined  experimentally  by  plotting  the  measured  values  of  q  and  Ap.  The  slope 
F  is  given  at  every  point  by: 

q 

F  =  — —  =  Tunnel  Correction  Factor  (5) 

Ap 

Once  F  has  been  determined,  test  section  velocity  can  be  calculated  from  Ap  by: 
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Calibration  measurements  were  taken  at  numerous  velocities  over  the  range  of  40  - 
140  ft  sec.  The  two  differential  pressures  were  recorded  at  each  calibration  point. 
Figure  10  shows  the  plot  of  q  versus  Ap  for  this  calibration.  The  relationship  is  ob¬ 
served  to  be  linear  and  the  tunnel  correction  factor  is  1.126.  This  value  differed  signif¬ 
icantly  from  the  calibration  factor  of  1.075  which  has  been  in  common  use  at  NPS  since 
1982.  The  author  speculates  that  the  original  1982  calibration  was  made  prior  to  in¬ 
stallation  of  the  reflection  plane.  Another  possibility  is  that  the  new  calibration  factor 
reflects  the  numerous  small  changes  made  to  the  tunnel  over  the  years. 

B.  HOT-WIRE  ANEMOMETER  CALIBRATION 

The  purpose  of  hot-wire  anemometer  calibration  was  to  determine  the  equation  re¬ 
lating  hot-wire  output  voltage  to  velocity.  This  equation  was  used  in  the  data  conver¬ 
sion  program  CONVDAT.FOR  described  earlier  in  Chapter  III.  The  hot-wire  was 
calibrated  before  and  after  each  data  run.  Hot-wire  calibration  consists  of  three  steps: 

1  FA- 100  system  setup  and  initialization 
Calibration  data  acquisition 
Curve  plotting  and  least  squares  fit 


Initialization  procedures  for  the  I  FA- 100  system  followed  the  steps  for  routine  op¬ 
eration  outlined  in  the  instruction  manual  {Ref.  8:  pp.  3-19J.  The  setup  parameters  for 
the  I  FA- 100  were  as  follows: 


Bridge 

High  Pass  Filter 
Low  Pass  Filter 
Voltage  Offset 
Voltage  Gain 


Standard  Bridge  Configuration 
Not  used 
20  KHz 
-1  volt 
5 


As  mentined  in  Chapter  II,  the  low  pass  filter  was  set  to  the  Nyquist  cutoff  frequency 
calculated  by  Equation  1.  The  Nyquist  cutoff  frequency  for  this  investigation  was  25 
KHz.  The  low  pass  filter  could  only  be  set  in  increments  of  10  KHz,  however,  so  the 
filter  was  set  to  20  KHz— the  closest  lower  frequency.  Voltage  offset  and  gain  were  used 
to  improve  the  measurement  resolution  of  the  analog-to-digital  converter  by  expanding 
the  output  voltage  range  of  the  hot-wire  anemometer,  as  discussed  in  Chapter  II. 

The  hot-wire  was  calibrated  against  a  reference  pitot-static  probe  located  near  the 
top  of  the  test  section.  The  hot-wire  was  traversed  to  a  position  adjacent  to  the  pitot- 
static  probe  for  each  calibration. 
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Figure  10.  Wind  Tunnel  Test  Section  Calibration 
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The  hot-wire  was  calibrated  at  or  near  seven  selected  velocities  (1,  2,  3,  4.  6.  8,  and 
10  centimeters  Ap  across  the  contraction  cone)  using  the  software  routine  HWCAL.BAS 
described  in  Chapter  II.  The  calibration  velocities  are  centered  about  the  expected  op¬ 
erating  speed  of  95  ft,  sec  and  extended  over  the  range  40  -  140  ft  sec. 

The  final  calibration  curve  was  obtained  by  plotting  hot-wire  voltage  squared  versus 
the  square  root  of  velocity.  A  nearly  straight  line  relationship  results  for  this  fourth  or¬ 
der  fit.  Most  graphing  routines  will  plot  the  data  and  provide  the  equation  of  the  linear 
least  squares  best  fit  line.  This  linear  equation  is  of  the  form: 

y  =  Ax  +  B  (7) 

where  y  =  v; velocity  and  x  =  (hot-wire  voltage)2.  Using  the  constants  A  and  B  from 
Equation  7,  the  hot-wire  calibration  equation  becomes: 

velocity  =  [A(voltage)2  +  B]2  (8) 

This  relationship  can  readily  be  used  in  any  data  conversion  program  such  as 
DATCOXV.FOR  (Chapter  II)  to  calculate  velocity. 

C.  HOT-WIRE  POSITIONING 

Great  attention  was  given  to  the  task  of  positioning  the  hot-wire  in  the  boundary 
layer  adjacent  to  the  wing.  The  objective  was  to  accurately  position  the  hot-wire  as  close 
as  possible  to  the  wing  surface  while  not  damaging  either  the  wing  or  the  hot-wire. 

Vaught  and  Fernholz  provide  an  excellent  summary  of  the  various  methods  in  use 
today  [Ref.  1 1].  The  optical  sighting  method  was  used  for  this  investigation.  Kindelspire 
also  used  this  technique  in  prior  work  at  NPS  [Ref.  4:  pp.  12-13]. 

Several  modifications  to  Kindelspire's  methods  were  made  by  this  investigator  to 
improve  positioning  accuracy  and  repeatability.  The  techniques  employed  are  described 
below. 

The  height  of  the  hot-wire  above  the  wing  surface  was  measured  by  observing  the 
hot-wire  and  its  reflection  in  the  wing  through  a  surveyor's  sighting  level.  The  sighting 
level  was  positioned  approximately  six  feet  behind  and  to  the  left  of  the  tunnel  test  sec¬ 
tion.  Figure  1 1  shows  the  sighting  level  in  position  adjacent  to  the  NPS  low  speed  wind 
tunnel. 

As  noted  previously  in  Chapter  II,  the  wing  was  finished  in  high  gloss  black  lacquer 
to  enhance  the  reflected  image  of  the  hot-wire  probe.  The  distance  between  the  hot-wire 
probe  and  its  image  was  determined  by  using  the  horizontal  cross-hair  and  vertical  offset 
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Figure  11.  Sun evor's  Sighting  Leu-l 

. ernier  or.  the  sighting  level.  I  lie  vortical  \ ernier  was  adjustable  cner  a  rar.ee  el  -  4' 
. err.ier  unit'  front  the  center  position.  One  \ ernier  unit  was  determined  to  be  1 1 2 
r.wi.e'  of  '.ertis.h  probe  di'p’.aeenvent  using  the  traverser,  with.  *'."<"•! 25  incites  stepping 
..o.;:\tp..  as  a  calibration  reference.  'I  he  sighting  level  could  be  repeated!;-  set  and  read 
.  1  2  •.  errucr  unit  result;:..-  m  a  positioning  accuracx  ol*  •  M  >*  *  1 1*  incite-. 

The  actual  probe  to  wing  distance  is  1  2  the  distance  observed  between  the  two  un¬ 
ices,  I  he  hot-wire  images  were  brought  to  within  three  vernier  units  separation  for  data 
.olieeted  at  the  chord  location  resulting  in  an  initial  hot-wire  distance  of  <U»<*4S  r 
'."i ms  inches.  I  he  probe  was  brought  to  two  vernier  units  separation  for  all  other  wing 
'o-ni;;>  .  'D  ..  and  eitordi  which  equated  to  an  initial  distance  of  2  r 


An  external  flood  ligiit  was  directed  at  the  wing  section  and  tite  hot-wire  probe 
urw.c  the  s;g;,t:nc  pro.es-  to  enhance  the  image  in  the  wtng.  1  he  polished  high  gloss 
2  a. 1.  on  the  w.nc  provided  a  clear  and  sharp  relleeted  imaee  cumparabic  q uni- 


Figure  12.  Hot-Wire  Probe  and  Reflection 

ity  to  the  actual  probe  as  shown  in  Figure  12.  No  correction  was  made  for  the  slight 
angle  of  the  wing  surface  due  to  camber  as  this  angle  was  less  than  five  degrees  at  the 
points  examined  resulting  in  a  positioning  error  of  less  than  three  percent. 

D.  ACQUISITION  OF  EXPERIMENTAL  DATA 

1.  Test  Emironnient  Conditions  and  Setup  Parameters 

All  test  were  conducted  at  a  chord  Reynolds  number  of  5(><UX)0.  Tunnel  ve¬ 
locity  for  a  given  Reynolds  number  is  given  by: 

R.u 

/ 1 1  i 


R.  =  Chord  Reynolds  Number 

l:  =  Viscosity 
/>  =  Density 
c  =  Chord  Length 


Solving  Equation  9  using  typical  NPS  ambient  conditions  (temperature  =  68  degrees  F, 
pressure  =  29.95  inches  Hg)  where  Rc  =  500.000  and  c=  10  inches  yields  a  tunnel  oper¬ 
ating  velocity  of  95  ft  sec.  All  data  runs  were  made  at  this  reference  free  stream  velocity. 

The  setup  parameters  that  remained  constant  throughout  the  investigation  are 
summarized  below: 


Wing  Angle  of  Attack  (AOA) 

0.71  degrees 

Wing  Chord 

10  inches 

Free  Stream  Reference  Velocity- 

95  ft  sec 

Chord  Reynolds  Number 

500,000 

Turbulence  Generator  Motor  Speed 

1,500  RPM 

Frequency  of  Turbulence  Pulses 

50  Hz 

Period  Between  Turbulence  Pulses 

20  milliseconds 

Free  Stream  Distance  between  Pulses 

22.8  inches 

Analog  Channels  Sampled 

2  channels 

Hot-Wire  Voltage 

Channel  1 

Timing  Data 

Channel  2 

Analog-to-Digital  Rate 

50  KHz  per  Channel 

Sample  Period 

0.92  seconds 

Number  of  Pulse  Periods  Sampled 

45 

Number  of  Samples  per  Pulse  period 

1000 

Period  Between  Each  Sample 

20  microseconds 

2.  Equipment  Configuration 

Figure  13  shows  the  equipment  configuration  used  for  this  investigation.  As 
noted  previously,  two  personal  computers  were  used  during  the  data  acquisition  phase. 
One  computer  was  dedicated  to  controlling  the  traverser  mechanism  while  the  other 
managed  data  collection  and  processing.  Two  computers  were  necessary  as  the  operat¬ 
ing  programs  and  data  (92,000  samples)  entirely  filled  the  free  memory  of  the  data  ac¬ 
quisition  computer. 

Data  collection  was  quite  routine  once  the  hot-wire  anemometer  was  set  up  and 
calibrated,  and  the  hot-wire  probe  initially  positioned  in  the  boundary  layer.  From  that 
point,  the  operator  had  only  to  monitor  the  operating  environment  (tunnel  velocity, 
tunnel  temperature,  turbulence  generator  motor  speed,  and  hot-wire  output),  step  the 
hot-wire  through  the  boundary  layer  using  one  computer,  and  collect  data  with  the  other 
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Figure  13.  Equipment  Configuration 

computer.  Approximately  seven  minutes  were  needed  to  collect  and  store  the  data  at 
any  one  position  in  the  boundary  layer.  While  the  actual  data  sampling  period  was  less 
than  one  second,  the  computer  took  over  5  minutes  to  transfer  the  data  from  random 
access  memory  to  permanent  disk  storage. 

3.  Flow  Characterization 

Data  were  collected  to  characterize  the  flow  in  the  following  regions  of  interest: 
Test  Section  Free  Stream 
Periodic  Turbulence  Pulses 

Boundary  Layer  Response  to  Periodic  Turbulence  Pulses 
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a.  Test  Section  Free  Stream 

Data  were  collected  in  the  free  stream  of  the  test  section  with  all  equipment 
removed  to  determine  the  intensity  and  spectral  content  of  the  flow.  This  step  served 
to  document  the  steady  state  baseline  characteristics  of  the  free  stream  flow  at  9C  ft  sec 
and  also  allowed  this  investigator  to  compare  the  free  stream  turbulence  data  obtained 
in  this  investigation  with  that  obtained  previously  by  other  investigators. 

b.  Periodic  Turbulence  Pulses 

Data  were  collected  with  just  the  periodic  turbulence  generator  installed  and 
operating  in  the  test  section.  The  tunnel  and  the  turbulence  generator  were  set  at  their 
normal  operating  parameters  of  95  ft  sec  and  50  Hz  respectively.  The  hot-wire  was  po¬ 
sitioned  at  the  approximate  vertical  position  of  the  upper  wing  surface  (approximately 
halfway  between  the  hub  and  the  rod  tip).  Measurements  were  taken  four,  six,  and 
eight  inches  downstream  from  the  rotating  rod.  These  three  positions  coincided  with  the 
30%,  50%,  and  70%  chord  locations  on  the  wing.  These  tests  served  to  document  the 
time-varying  velocity,  turbulence,  and  spectral  characteristics  of  the  turbulence  pulse 
generated  by  the  rotating  bar. 

c.  Boundary  Layer  Response  to  Periodic  Turbulence  Pulses 

Measurements  were  made  in  the  boundary  layer  at  the  30%,  50%.  and  70% 

chord  locations  on  the  upper  surface  of  the  wing  with  the  periodic  turbulence  generator 
operating  at  50  Hz.  The  wing  AOA  was  measured  at  0.71  degrees  and  the  free  stream 
velocity  was  set  at  95  ft  sec.  Data  were  taken  at  up  to  50  points  over  the  estimated 
boundary  layer  thickness  at  each  chord  location.  The  spacing  of  the  sample  points  taken 
at  each  of  the  three  chord  locations  is  given  in  Appendix  B. 

All  data  measurements  were  taken  over  a  four  day  period.  One  additional 
day  of  tunnel  operation  was  needed  following  the  data  collection  runs  in  order  to  inves¬ 
tigate  the  hot-wire  drift  problems  addressed  in  Chapter  IV. 
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IV.  RESULTS 


A.  HOT-WIRE  ANEMOMETER  CALIBRATION 

Hot-wire  calibration  was  conducted  before  and  after  each  boundary  layer  investi¬ 
gation?  to  confirm  the  stability  of  the  hot-wire  anemometry  system  and  to  ensure  the 
accuracy  of  the  data.  Roane  [Ref.  12:  pp.  40-41]  experienced  some  calibration  curve  drift 
in  previous  hot-wire  work  in  the  NPS  low  speed  wind  tunnel.  As  such,  a  minor  amount 
of  drift  was  expected.  This  investigator  planned  to  average  the  initial  and  final  cali¬ 
bration  curves  and  use  the  resulting  average  to  convert  each  day's  data.  The  tabulated 
data  from  the  calibrations  runs  on  each  day  tended  to  confirm  the  expectation  of  only 
slight  drift  and.  therefore,  the  data  acquisition  runs  were  completed  as  planned.  Subse¬ 
quent  detailed  analysis  of  the  calibration  data  for  day  three,  however,  showed  that  the 
initial  and  final  curves  differed  by  approximately  14%.  The  drift  for  day  two  was  less 
but  still  significant  at  approximately  8%.  The  drift  for  day  one  could  not  be  determined 
as  the  final  calibration  data  were  lost.  The  results  of  the  calibration  runs  for  day  three, 
presented  in  Figure  14,  show  the  divergence  between  the  initial  and  final  calibration 
data.  The  magnitude  of  the  error  was  too  great  to  resolve  by  averaging.  In  addition, 
the  cause  of  the  drift  and  its  characteristics  were  unknown. 

This  investigator  planned  and  conducted  an  additional  six-hour  tunnel  run  on  the 
following  day  in  order  to  document  and  characterize  the  observed  hot-wire  drift.  The 
tunnel  was  run  at  95  ft  sec  with  just  the  pitot-static  and  hot-wire  probes  installed  in  the 
test  section.  Calibrations  were  conducted  at  the  start  and  every  hour  thereafter.  The 
results  of  this  run.  plotted  in  Figure  15.  show  the  shift  of  the  calibration  curve  over  time. 
The  drift  rate  is  greatest  at  the  beginning  and  becomes  very  small  after  six  hours.  The 
effects  of  the  drift  and  the  induced  error  are  greatest  at  the  higher  velocities  due  to  the 
fourth  order  relationship  that  exists  between  hot-wire  voltage  and  velocity. 

The  drift  rate  was  graphically  analyzed  at  a  reference  velocity  of  100  ft,  sec  and  ex¬ 
panded  by  extrapolation  to  cover  an  eight-hour  period.  Figure  16  shows  the  change  in 
calibration  of  the  hot-wire  over  eight  hours.  Notice  that  nearly  50%  of  the  drift  occurs 
in  the  first  two  hours  of  operation, 

5  The  final  calibration  data  for  day  one  was  not  saved  due  to  operator  error. 
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Figure  16.  Hot-Wire  Calibration  Drift  Over  Eight  Hour  Period 
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The  autnor  applied  the  drift  vs.  time  relationship  of  Figure  16  to  each  day's  cali¬ 
bration  data  in  order  to  generate  a  set  of  time  varying  calibration  equations.  The 
equations  were  staggered  such  that  the  maximum  velocity  difference  between  sequential 
calibration  curves  was  one  ft  sec  or  less.  These  corrections  assumed  that  the  drift  rate 
determined  during  the  test  run  was  representative  of  the  drift  experienced  on  each  day 
and  that  the  run  duration  for  each  days  was  eight  hours.  Corrections  were  applied  only 
to  the  data  taken  on  day  two  (50%  chord)  and  day  three  (30%  chord). 

This  investigator  believes  the  hot-wire  drift  is  a  result  of  environmental  contam¬ 
ination  of  the  hot-wire  probe.  Hot-wire  contamination  appears  to  be  a  well  known 
phenomenon  among  experienced  researchers,  but  a  subject  about  which  little  is  printed 
in  current  hot-wire  operating  manuals  or  in  fluid  flow  measurement  texts.  An  older 
publication  by  TSI  has  this  to  say  about  hot-wire  stability: 

For  a  mean  flow  measurement  to  be  useful  it  must  maintain  its  calibration  over 
long  time  periods.  Hot  wires  are  notoriously  troublesome  in  this  regard  due  to  their 
small  size  making  them  susceptible  to  picking  up  small  contaminants  ....  The 
larger  hot  film  sensor  has  been  a  great  improvement  in  this  regard.  The  effects  of 
surface  contamination  tend  to  be  a  function  of  sensor  size  and  a  ceramic  coating  on 
the  sensor  surface  prevents  corrosion  and  erosion.  Hence,  the  hot  film  is  far  supe¬ 
rior  for  mean  measurements.  [Ref.  13:  p.  91] 

Figure  17.  reprinted  from  the  above  TSI  publication,  shows  the  typical  change  in 
calibration  of  a  hot-wire  sensor  and  a  hot-film  sensor  vs.  time  when  subjected  to  an  ar¬ 
tificially  dirty  flow  stream.  The  drift  vs.  time  curve  for  this  investigation  (Figure  16)  is 
remarkably  similar  to  the  TSI  drift  curve.  With  this  knowledge,  it  was  not  difficult  to 
pinpoint  the  source  of  contamination  to  the  wind  tunnel  itself. 

For  the  month  prior  to  this  investigation,  the  tunnel  was  used  for  a  series  of  flow 
visualization  studies  using  an  oil-based  smoke  as  a  visualization  medium.  This  type  of 
medium  tends  to  produce  a  certain  amount  of  contaminants  which,  due  to  the  closed 
circuit  design  of  the  tunnel,  were  trapped  in  the  tunnel  environment.  This  oil  residue  is 
the  likely  contaminant  and  cause  of  the  observed  drift. 

Little  can  be  done  about  the  tunnel  contamination,  but  one  likely  solution  to  the 
drift  problem  is  the  use  of  a  hot-film  sensor  instead  of  a  hot-wire  probe.  As  noted  pre¬ 
viously,  the  hot-film  is  virtually  impervious  to  particulate  contamination.  A  second 
solution  would  be  to  pre-run  the  tunnel  prior  to  taking  data  for  six  to  eight  hours  to 
allow  the  hot-wire  to  stabilize. 
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Figure  17.  Change  in  Calibration  of  Hot-Wire  vs.  Hot- Film  When  Subjected  to  a 
Dirty  Flow  Stream:  [Ref.  13:  p.  92] 

B.  ANALYSIS  OF  TIME  VARYING  DATA 

Ensemble  averaging  methods  were  employed  to  characterize  the  non-stationary,  pe¬ 
riodic  data  acquired  during  this  investigation.  The  periodic  nature  of  the  data  made  it 
possible  to  collect  a  series  of  repetitive  time  histories,  or  pulse  passage  cycles,  each  as¬ 
sociated  with  the  passage  of  one  turbulence  pulse.  This  series  of  time  histories  is  called 
an  ensemble.  A  typical  ensemble  consisted  of  a  series  of  45  pulse  passage  cycles.  The 
number  of  passage  cycles  in  one  ensemble  was  essentially  limited  by  available  computer 
memory.  For  this  investigation,  45  cycles  was  the  upper  limit.  Howard  [Ref.  6:  pp. 
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114-115]  showed  that  20  or  more  pulses  were  needed  in  each  ensemble  to  accurately 
portray  the  actual  mean  values,  substantiating  that  45  pulses  were  sufficient  in  this  case. 

1.  Ensemble  Characteristics 

Each  of  the  45  pulse  passage  cycles  making  up  the  typical  ensemble  for  any  one 
data  point  is  slightly  different  due  to  the  random  nature  of  the  turbulent  components 
of  the  signal.  Figure  18  and  Figure  19  show  two  sequential  pulse  passage  cycles  {in¬ 
stantaneous  velocity  vs.  sample  number)  from  a  typical  ensemble  {Test  Case  A)  taken 
near  the  wing  surface  at  70%  chord.  These  data  were  collected  using  a  single  hot-wire 
sensor.  The  velocity  data  were  sampled,  or  digitized,  at  50  KHz  with  the  turbulence 
generator  operating  at  50  Hz.  This  digitizing  rate  results  in  approximately  1000  data 
points  per  passage  cycle.  The  period  of  each  passage  cycle  being  20  milliseconds,  each 
sample  represents  20  microseconds  in  time.  Time  is  not  recorded  explicitlv  but  can  be 
calculated  from  the  digitizing  rate  and  sample  number.  The  data  are  plotted  with 
straight  line  segments  between  each  of  the  sample  points. 

Each  passage  cycle  is  typically  made  up  of  three  representative  flow  regimes  as 
depicted  in  Figure  IS: 

Undisturbed  Flow 
Turbulence  Pulse 
Recovery  Period 

The  specific  characteristics  of  these  regimes  are  very  dependent  on  chord  location  and 
position  in  the  boundary  layer;  but  the  following  general  characteristics  hold  true  for 
each.  The  undisturbed  flow  segment  is  representative  of  the  natural  steady-state  flow  for 
the  wing  when  operated  in  an  undisturbed  free  stream.  Howard  [Ref.  6:  p.  81]  confirmed 
this  relationship  experimentally  by  comparing  the  drag  of  the  undisturbed  segment  with 
that  for  natural  flow  and  finding  the  two  equal.  The  turbulence  pulse  is  clearly  identified 
as  a  highly  turbulent  velocity  deficit.  The  recovery  period  is  broadly  defined  as  the  pe¬ 
riod  from  the  end  of  the  turbulence  pulse  to  the  point  where  the  flow  returns  to  the  un¬ 
disturbed  state. 

2.  Mean  Velocity  Determination 

The  author  employed  program  MEANSMU.FOR  described  in  Chapter  II  to 
extract  the  slowly  varying  mean  velocity  from  the  ensemble  of  instantaneous  velocities. 
In  principle,  MEANSMU  filters  out  the  high  frequency  turbulence  components  while 
leaving  the  low  frequency  variations  of  the  signal  mean  intact.  The  amount  of  smooth- 
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Figure  18.  Pulse  Cycle  Number  One  for  Test  Case  A:  (Fs  =  50  KHz) 
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Figure  19.  Pulse  Cycle  Number  Tuo  for  Test  Case  A:  (Fs  =  50  KHz) 
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ing  is  determined  by  the  user  selectable  cutofT  frequency  of  the  digital  low  pass  filter  in 
program  MEANS  Ml'. 

Program  MEANSMU  extracts  the  mean  through  a  two  step  process.  First,  the 
45  time  histories  are  ensemble  averaged  to  generate  a  single  "representative"  time  his¬ 
tory.  Ensemble  averaging  consists  of  simply  stacking  the  45  time  histories  one  on  an¬ 
other  and  taking  the  average  of  the  45  values  at  each  sample  point  in  the  time  history. 
Figure  20  presents  the  ensemble  average  of  the  45  pulse  passage  cycles  for  Test  Case 
A.  Note  that  the  large  turbulent  fluctuations  observed  in  Figure  18  and  Figure  19  have 
been  reduced  but  not  entirely  eliminated.  The  three  flow  regimes  are  also  quite  pro¬ 
nounced  and  easily  identifiable. 

Figure  21  through  Figure  24  show  the  application  of  digital  low  pass  filtering 
to  the  ensemble  average  of  Figure  20  with  low  pass  filter  cutoff  frequencies  (fc)  of  1082 
Hz,  541  Hz.  360  Hz,  and  267  Hz  respectively.  Note  that  the  amount  of  smoothing  is 
directly  related  to  fc ;  the  lower  fc  the  greater  the  smoothing.  360  Hz  was  choosen  as  the 
optimal  cutoff  frequency  based  on  the  findings  of  Landrum  and  Macha  [Ref.  14]  in  their 
analysis  of  the  spectra  documenting  the  growth  of  turbulent  disturbances  through  the 
transition  process.  They  found  that  the  free  stream  energy  is  contained  below  200  Hz 
while  the  turbulence  energy  is  generally  contained  above  550  Hz  in  the  early  transition 
stages.  The  cutoff  frequency  was  choosen  to  lie  between  these  two  frequencies. 
Coincidentially,  the  fc=  360Hz  curve  is  relatively  smooth  and  fiat  throughout  the  undis¬ 
turbed  flow  regime  and  displays  good  peak  response  over  both  the  turbulence  pulse  and 
recovery  period  regimes. 

3.  Turbulence  Level  Determination 

Turbulence  values  of  the  time- varying,  periodic  signals  were  calculated  using  the 
program  TURBIN. FOR  described  in  Chapter  II.  TURBIN  accomplishes  this  by  calcu¬ 
lating  the  instantaneous  turbulence  at  each  point  in  each  of  the  45  time  histories  and 
then  ensemble  averaging  the  instantaneous  turbulence  magnitudes.  The  result  is  a 
time-varying  representation  of  the  mean  turbulence  at  each  sample  point  in  the  pulse 
period.  It  is  important  to  note  that  TURBIN  also  employs  user  selectable  digital  low 
pass  filtering  in  its  routines.  The  author  used  360  Hz  in  this  routine  to  maintain  conti¬ 
nuity  among  the  various  output  data. 

Figure  25  shows  the  turbulence  trace  for  Test  Case  A.  Note  that  the  represen¬ 
tative  flow  regimes  (undisturbed  flow,  turbulence  pulse,  and  recovery  period)  are  as 
clearly  identifiable  here  as  in  the  mean  velocity  traces  in  the  previous  section. 
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Figure  20.  Ensemble  Average  for  Test  Case  A:  (Fs=  50  KHz) 
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Figure  21.  Smoothed  Ensemble  Average  for  Test  Case  A,  Fc=  1082  Hz:  (Fs=  50 
KHz) 
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Figure  22.  Smoothed  Ensemble  Average  for  Test  Case  A,  Fc=54I  Hz:  (Fs=50 
KHz) 


41 


100  200  300  400  500  600  700  800  900  1000 

Sample  Number 


100 


(oas/u)  A^popA 


Figure  23.  Smoothed  Ensemble  Average  for  Test  Case  A,  Fc=360  Hz:  (Fs  =  50 
KHz) 
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Figure  24.  Smoothed  Ensemble  Average  for  Test  Case  A,  Fc=267  Hz:  (Fs  =  50 
KHz) 
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Figure  25.  Smoothed  Mean  Turbulence  for  Test  Case  A:  (Fc=  360  Hz,  Fs  =  50 
KHz) 
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4.  Spectral  Analysis 

Spectral  analysis  of  the  time-varying,  periodic  signals  was  conducted  using  pro¬ 
gram  ENSPSDAV.  FOR  described  in  Chapter  11.  ENSPSDAV  performs  spectral  analy¬ 
sis  on  a  user  selectable  256  sample-wide  window  of  the  pulse  passage  cycle  time  history. 
In  theory,  spectral  analysis  can  only  be  accomplished  on  a  discrete  (digitized)  time  series 
of  infinite  length.  In  practice,  the  window  width,  or  number  of  samples,  can  be  reduced 
through  the  application  of  special  techniques  to  the  Fast  Fourier  Transform-the  com¬ 
promise  being  loss  of  frequency  resolution  with  reduction  in  window  width.  The  256 
sample-wide  window  used  in  this  investigation  was  selected  as  the  best  balance  of  these 
factors.  Each  flow  regime  fits  into  the  window  and  the  resulting  frequency  resolution  is 
190  Hz.  The  reader  must  keep  in  mind  that  the  output  data  from  ENSPSDAV  repres¬ 
ents  the  average  taken  over  the  entire  window.  Thus,  dynamic  changes  that  occur  in  the 
How  regime  over  the  window  period  will  be  lost  in  the  resulting  average. 

Figure  26  and  Figure  27  show  the  output  from  program  ENSPSDAV  for  the 
turbulence  pulse  regime  of  Test  Case  A.  Figure  26  displays  the  velocity  spectrum  of  the 
signal.  The  vertical  scale  is  a  dB  measurement  of  the  relative  turbulence  energy  with  0 
dB  representing  the  highest  magnitude  component  of  the  signal.  The  horizontal  scale  is 
frequency  in  Hz  plotted  out  to  the  Nyquist  frequency.  Each  10  dB  reduction  in  the  ve¬ 
locity  spectrum  represents  a  10  fold  reduction  in  power  level.  For  example,  the  power 
level  at  4.500  Hz  (-20  dB  )  is  1  100  the  level  of  the  highest  magnitude  component.  The 
mean  has  been  removed  from  the  signal^  so  that  only  the  spectrum  above  0  Hz  is  ana¬ 
lyzed.  It  is  important  to  note  that  the  spectrum  display  provides  no  information  about 
the  absolute  power  contained  in  the  spectrum,  but  only  displays  power  levels  relative  to 
the  highest  magnitude  in  the  spectrum.  Figure  27  displays  the  total  (cumulative)  power 
distribution  plotted  with  regard  to  frequency.  The  total  power  curve  asymptotically 
approaches  the  value  of  the  total  power  in  the  spectrum.  The  shape  of  this  curve  is  in¬ 
dicative  of  the  power  distribution  within  the  spectrum.  The  sharp  knee  observed  in  this 
plot  indicates  that  the  low  frequency  components  of  the  signal  contain  most  of  the 
spectral  power  with  very  little  power  in  the  higher  frequency  components. 

C.  FREE  STREAM  ANALYSIS 

The  results  of  the  free  stream  analysis  are  presented  for  flow  at  95  ft  sec  with  all 
equipment  removed  from  the  test  section.  Figure  28  is  a  plot  of  velocity  response  over 

6  The  mean  is  removed  for  all  the  velocity  spectrum  and  total  power  distribution  plots  in  this 
report. 
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Figure  26.  Velocity  Spectrum  for  Test  Case  A,  Turbulence  Pulse 

1000  samples.  The  signal  is  flat  with  a  high  frequency  turbulent  component  of  small 
amplitude  superimposed  on  the  mean.  The  free  stream  turbulence  was  calculated  at 
0.182%  by  the  data  analysis  methods  of  this  investigation.  This  figure  compares  favor¬ 
ably  with  the  free  stream  turbulence  intensity  of  0.23%  observed  by  Roane  (Ref.  12:  p. 
41]  at  a  velocity  of  116.4  ft.'sec.  The  slightly  lower  value  of  free  stream  turbulence  ob¬ 
served  here  is  not  unexpected  due  to  the  slower  free  stream  velocity. 

Figure  29  shows  the  velocity  spectrum  of  the  free  stream.  The  curve  shows  a  sharp 
drop  in  intensity  with  increasing  frequency  indicating  that  the  major  frequency  compo¬ 
nents  are  at  the  lower  frequencies. 
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Figure  27.  Total  Power  Distribution  for  Test  Case  A,  Turbulence  Pulse 

The  total  power  distribution  of  the  Tree  stream  is  depicted  in  Figure  30.  The  rapid 
rise  and  sharp  knee  of  this  curve  confirm  the  concentration  of  power  in  the  lower  fre¬ 
quencies.  The  frequency  components  above  2  KHz  provide  little  additional  power  as 
shown  by  the  flatness  of  the  curve  above  this  point. 

D.  TURBULENCE  PULSE  ANALYSIS 

The  turbulence  pulse  was  examined  at  three  positions  (4",  6",  and  8")  downstream 
of  the  spinning  rod  with  the  wing  and  wingmount  assembly  removed.  The  smoothed 
mean  velocity  time  histories  for  the  three  positions  are  shown  in  Figure  31.  The  plots 
show  a  nearly  symmetric  velocity  deficit  trace  at  all  three  positions.  Pulse  duration  is 
roughly  250  samples  or  5  milliseconds  and  the  deficit  maxima  are  constant  at  -20  ft,  sec. 
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Figure  28.  Test  Section  Free  Stream  Flow:  (Fs  =  50KHz) 
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Figure  29.  Free  Stream  Velocity  Spectrum 

1  he  width  of  the  deficit  appears  to  increase  slightly  with  increasing  distance  downstream. 
The  velocity  deficit  trace  observed  for  the  spinning  rod  is  very  similar  to  that  generated 
by  a  stationary  cylinder  in  steady  cross  flow.  This  similarity  shows  that  spinning  the  rod 
docs  not  significantly  alter  the  classical  response  of  the  cylinder. 

Figure  32  is  the  plot  of  the  smoothed  mean  turbulence  traces  of  the  turbulence  pulse 
at  the  same  positions.  The  pulse  is  intensely  turbulent  with  a  peak  value  of  10.5%.  The 
turbulence  traces  coincide  with  the  velocity  traces;  however,  the  turbulence  traces  be¬ 
come  somewhat  non-symmctric  downstream.  In  addition,  the  turbulence  appears  to 
anticipate  or  lead  the  velocity  deficit  by  up  to  60  samples  or  1.2  milliseconds.  Most 
notably,  the  turbulence  magnitude  decreases  quite  rapidly  with  distance  downstream. 
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Figure  30.  Free  Stream  Total  Power  Distribution 

This  decay  is  not  unexpected  and  is  generally  attributed  to  dissipation  of  the  high  fre¬ 
quency  turbulence  components  due  to  viscosity  [Ref.  15:  p.  48). 

Figure  33  shows  the  velocity  spectra  of  the  turbulence  pulse  regime  and  the  undis¬ 
turbed  flow  regime  taken  at  4"  downstream  from  the  spinning  rod.  The  turbulence  pulse 
trace  decays  more  rapidly  and  to  a  lower  level  than  the  undisturbed  flow  trace.  Addi¬ 
tionally,  the  undisturbed  flow  trace  levels  out  at  approximately  -20  dB  for  frequencies 
above  5  KHz  while  the  turbulence  pulse  trace  bottoms  out  at  -42  dB.  These  plots  show 
the  turbulence  pulse  to  be  a  low  frequency  dominated  event  while  the  undisturbed  flow 
regime  has  components  at  much  higher  frequencies. 
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Figure  31.  Turbulence  Pulse  Mean  Velocity:  (Fs=50KHz) 
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Figure  32.  Turbulence  Pulse  Mean  Turbulence:  (Fs  =  50KHz) 
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Figure  33.  Turbulence  Pulse  Velocity  Spectrum 

The  total  power  curves  for  the  two  regimes  are  plotted  in  Figure  34.  The  intensity 
of  the  turbulence  pulse  is  clearly  shown  here:  the  pulse  power  is  4  orders  of  magnitude 
greater  than  the  undisturbed  How.  The  rapid  rise  and  sharp  knee  of  the  pulse  trace  also 
confirms  the  intense,  low  frequency  nature  of  the  turbulence  pulse.  Very  little  power  is 
added  by  the  frequency  components  above  1  KHz. 

E.  BOUNDARY  LAYER  ANALYSIS 

The  response  of  the  boundary  layer  to  periodic  turbulence  was  investigated  at  three 
locations  on  the  upper  surface  of  the  wing:  30%  chord,  50%  chord,  and  70%  chord. 
These  three  locations  were  chosen  as  the  natural  steady-state  (undisturbed)  flow  was 
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Figure  34.  Turbulence  Pulse  Total  Poner  Distribution 


54 


laminar  at  30%  chord,  transitional  turbulent  at  50%  chord,  and  fully  turbulent  at  70'% 
chord. 

1.  30%  Chord 

The  ensemble-average  velocity  and  turbulence  profiles  are  shown  in  Figure  35. 
These  plots  consist  of  a  series  of  time-sequenced  profiles,  or  "snapshots",  taken  at  several 
points  in  the  time  history.  These  profiles  depict  the  dynamic  characteristics  of  the 
boundary  layer  throughout  one  pulse  passage  cycle. 

The  profiles  show  that  the  flow  at  30%  chord  transitions  from  undisturbed 
laminar  flow  with  very  little  turbulence  (sample  100);  to  highly  turbulent  flow  with  a 
nearly  constant  level  of  turbulence  at  the  turbulence  pulse  center  (sample  275);  and  back 
again  to  laminar  flow  (sample  550).  It  is  interesting  to  note  how  quickly  and  uniformly 
the  boundary  layer  responds  to  the  turbulence  pulse.  The  turbulence  profiles  show  a 
nearly  uniform  turbulence  increase  and  corresponding  decrease  with  the  passing  of  the 
turbulence  pulse.  At  the  same  time,  the  velocity  profiles  transition  quickly  from  laminar 
to  turbulent  and  back  to  laminar  in  step  with  the  turbulence  pulse.  The  velocity  increase 
observed  near  the  surface  of  profile  350  is  unusual.  The  velocity  increase  occurs  just 
after  the  turbulence  pulse  and  is  still  evident  in  the  following  profile  (sample  400).  The 
possible  reasons  for  this  behavior  will  be  examined  later  in  this  section. 

The  smoothed  mean  velocity  time  histories  and  the  smoothed  mean  turbulence 
time  histories  of  the  pulse  passage  cycle  are  shown  in  Figure  36  and  Figure  37  respec¬ 
tively.  These  traces  represent  the  time  histories  at  increasing  vertical  distances  (y)  above 
the  wing  surface.  The  bottom  trace  shows  the  flow  history  taken  closest  to  the  wing 
surface  (near-surface  region)  while  the  top  trace  depicts  external  flow  outside  of  the 
boundary  layer. 

The  velocity  and  turbulence  traces  in  the  external  flow  region  are  typical  of  the 
response  observed  downstream  of  the  periodic  turbulence  generator  (Figure  31  and 
Figure  32);  i.e.,  a  symmetric,  bell-shaped  pulse  followed  by  an  immediate  return  to  the 
free  stream  or  undisturbed  state. 

The  velocity  time  history  traces  show  a  velocity  deficit  of  1 1  ft  sec  in  the  ex¬ 
ternal  flow  due  to  the  turbulence  pulse.  This  deficit  increases  as  the  surface  is  ap¬ 
proached  until  very  close  to  the  wing.  Quite  near  the  surface,  the  velocity  deficit 
decreases  significantly  and  is  followed  by  a  sharp  velocity  increase  just  after  the  turbu¬ 
lence  pulse.  This  velocity  increase  remains  for  some  time  and  is  slow  to  dissipate.  The 
turbulence  traces,  on  the  other  hand,  show  that  the  turbulence  intensity  of  the  pulse  is 
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Figure  37.  Mean  Turbulence  Time  Histories  in  the  Boundary  Layer,  30% 
chord:  (Ts  =  50KHz) 
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felt  uniformly  throughout  the  boundary  layer  with  only  a  slight  reduction  in  magnitude 
and  stretching  of  the  pulse  very  near  the  surface. 

Figure  38  through  Figure  42  show  the  velocity  spectrum  and  total  power  dis¬ 
tribution  plots  for  the  three  major  flow  regimes  (undisturbed  flow,  turbulence  pulse,  and 
recovery  period)  at  points  in  the  boundary  layer.  Recall  that  the  velocity  spectrum  plot 
display  is  relative  to  the  highest  component  in  the  spectrum  while  the  total  power  dis¬ 
tribution  plot  depicts  the  cumulative  actual  power  with  regard  to  frequency.  Recall  also 
that  the  mean  has  been  removed  from  the  spectral  analysis.  Hence,  power  changes  in 
the  mean  due  to  mean  velocity  increases  or  deficits  are  not  reflected  in  the  spectral 
curves.  One  must  view  the  mean  velocity  and  turbulence  time  histories  along  with  the 
profiles  and  the  spectral  plots  as  a  group  in  order  to  understand  the  boundary  layer  be¬ 
havior. 

Notice  how  little  the  turbulence  pulse  affects  the  spectral  distribution  very  near 
the  surface  in  Figure  38.  The  velocity  spectra  for  the  undisturbed  flow  and  recovery 
period  regimes  are  very  similar  and  the  total  power  distribution  curves  for  the  two  re¬ 
gimes  are  comparable  in  shape  and  magnitude.  In  addition,  there  is  no  unusual  recovery 
period  in  this  region.  The  turbulence  pulse  comes  and  goes  without  affecting  the  spectral 
characteristics  of  the  flow  very  near  the  surface.  Higher  in  the  boundary  layer, 
Figure  40  and  Figure  41,  the  high  frequency  components  become  more  significant  in 
the  undisturbed  flow  and  recovery  period  regimes  while  the  total  power  contained  in 
both  regimes  drops  sharply.  This  indicates  more  spectral  power  (turbulence)  is  con¬ 
tained  in  the  boundary  layer  very  near  the  surface  than  in  the  intermediate  regions  or 
even  in  the  external  flow  region.  Note:  The  author  is  unable  to  explain  the  unusual 
shapes  of  the  undisturbed  flow  velocity  spectra  in  Figure  41  and  Figure  42  or  the  ap¬ 
parent  increase  in  undisturbed  flow  total  power  in  Figure  42. 

Together,  the  velocity  spectra  and  total  power  distribution  plots  show  that  the 
flow  recovers  very  quickly  from  the  turbulence  pulse  and  the  pulse  does  not  significantly 
affect  the  characteristics  of  the  recovery  period.  In  other  words,  the  undisturbed  flow 
and  recovery  period  velocity  spectra  and  total  power  distribution  curves  are  virtually 
identical  at  all  positions  in  the  boundary  layer. 

The  author  believes  an  understanding  of  the  mechanics  operating  in  the 
boundary  layer  at  this  chord  location  is  essential  in  understanding  the  more  complex 
responses  observed  at  50%  and  70%  chord.  The  discussion  that  follows  is  the  author's 
interpretation  of  the  operating  mechanics. 
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Figure  38.  Velocity  Spectra  and  Total  Potter  Distribution,  30%  Chord,  y  =  0.0032" 
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The  data  show  a  normally  laminar  flow  being  uniformly  disturbed  by  a  strong, 
highly  turbulent  pulse  from  the  external  flow.  The  boundary  layer  responds  to  the 
passing  pulse  by  momentarily  transitioning  to  extremely  turbulent  flow.  Once  the  pulse 
has  passed,  the  boundary  layer  rapidly  returns  to  low  turbulence  laminar  flow.  Note 
however,  that  the  velocity  time  history  graph  (Figure  36)  shows  that  the  mean  velocity 
quite  near  the  surface  increases  in  response  to  the  turbulence  pulse.  Furthermore,  this 
velocity  increase  remains  in  evidence  for  a  considerable  period  after  the  velocity  profile 
returns  to  laminar.  The  cause  of  this  velocity  increase  is  important. 

The  turbulence  time  histories  show  there  is  considerable  turbulence  and  mixing 
throughout  the  entire  boundary  layer  in  response  to  the  pulse.  This  seems  to  indicate 
that  momentum  transfer  down  through  the  bound  -.y  layer  is  possibly  responsible  for  the 
velocity  increase  near  the  surface.  The  velocity  time  histories  seem  to  support  this 
supposition  with  regard  to  observed  characteristics  of  the  velocity  deficit.  The  growth 
of  the  velocity  deficit  as  the  surface  is  approached  is  consistent  with  momentum  transfer 
down  through  the  layer. 

The  boundary  layer  region  very  near  the  wing  surface  cannot  transfer  momen¬ 
tum  due  to  the  solid  wing  boundary.  Thus,  the  near-surface  region  absorbs  all  the 
transferred  momentum  as  a  velocity  increase.  This  is  evidenced  by  the  sharp  decrease 
of  the  velocity  deficit  and  rapid  velocity  increase,  or  flow  acceleration,  just  following  the 
pulse.  These  two  phenomena  are  clearly  observed  in  the  bottom  two  velocity  time  his¬ 
tories.  Finally,  the  velocity  increase  remains  in  the  near-surface  region  due  to  the  posi¬ 
tive  velocity  gradient  ( )  in  the  boundary  layer  and  is  slowly  dissipated  by  the 
viscous  effects  of  the  laminar  layer. 

In  summary,  the  velocity  increase  in  the  near-surface  region  of  the  boundary 
layer  appears  to  be  the  result  of  a  downward  transfer  of  momentum  within  the  boundary 
layer  induced  by  the  high  turbulence  of  the  turbulence  pulse. 

2.  50%  Chord 

Figure  43  shows  the  ensemble  average  velocity  and  turbulence  profiles  at  50% 
chord.  The  boundary  layer  response  at  this  position  is  quite  unlike  the  simple  laminar- 
turbulent-laminar  response  observed  at  30%  chord.  To  begin  with,  the  undisturbed  flow 
at  sample  150  is  transitional  turbulent  with  5%  to  6%  turbulence  near  the  surface.  The 
flow  becomes  extremely  turbulent  at  sample  350  with  a  high  level  of  turbulence  across 
the  entire  profile.  During  the  recovery  period,  the  flow  becomes  what  appears  to  be  fully 
laminar  by  sample  650;  the  turbulence  is  extremely  low  over  the  entire  region  and  the 
boundary  layer  thickness  has  decreased  from  around  0.08  inches  at  sample  150  to  0.04 
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inches.  Samples  750  and  S50  show  the  flow  returning  to  the  original 
transitional  turbulent  state  as  evidenced  by  increasing  levels  of  turbulence  in  the  near 
surface  region. 

The  smoothed  velocity  time  histories  (Figure  44)  clearly  show  the  influence  of 
the  velocity  deficit  over  the  entire  boundary  layer.  The  corresponding  smoothed  turbu¬ 
lence  time  histories  (Figure  45),  however,  do  not  show  any  turbulence  penetration  into 
the  near-surface  region.  Apparently  the  high  level  of  natural  near-surface  turbulence 
overshadows  the  effects  of  the  turbulence  pulse.  The  apparent  answer  to  this  dilemma 
lies  in  the  velocity  spectra  and  total  power  distribution  plots  of  the  undisturbed  flow  re¬ 
gime  and  the  turbulence  pulse  regime  near  the  surface. 

Shown  in  Figure  46  are  the  velocity  spectra  and  total  power  distribution  curves 
taken  near  the  bottom  of  the  boundary  layer.  This  figure  corresponds  to  the  near¬ 
surface  region  of  the  mean  velocity  and  mean  turbulence  traces.  The  undisturbed  flow 
total  power  trace  is  quite  broad  banded  with  significant  power  contributions  from  the 
frequency  components  up  to  4  KHz  due  to  the  high  levels  of  turbulence  generated  near 
the  surface.  The  turbulence  pulse,  on  the  other  hand,  is  a  low  frequency  dominated 
phenomenon  having  no  significant  spectral  power  components  above  1  KHz.  Thus, 
while  the  turbulence  intensity  of  the  undisturbed  flow  and  turbulence  pulse  are  almost 
equal  in  magnitude,  the  spectral  makeup  of  each  regime  varies  greatly.  The  distinct  na¬ 
ture  of  the  turbulence  pulse  makes  it  felt  throughout  the  turbulent  near-surface  layers. 
The  effect  of  the  turbulence  pulse  on  the  undisturbed  flow  is  clearly  seen  in  this  figure. 
The  undisturbed  flow  spectrum  is  altered  by  the  turbulence  pulse  such  that  the  recovery 
period  spectrum  resembles  the  turbulent  pulse  spectrum. 

Figure  47  through  Figure  51  depict  the  changes  in  the  velocity  spectra  and  total 
power  distribution  at  increasing  vertical  distances  in  the  boundary  layer.  Near  the  sur¬ 
face  (Figure  46),  the  turbulence  pulse  clearly  alters  the  recovery  period  spectrum  by  re¬ 
moving  the  high  frequency  components  contained  in  the  undisturbed  flow  regime.  In 
the  middle  of  the  boundary  layer  (Figure  49),  the  recovery  period  spectrum  begins  to 
resemble  that  of  the  undisturbed  flow  showing  the  decreasing  impact  of  the  turbulence 
pulse  on  the  flow.  Finally,  in  the  external-flow  region  (Figure  51).  the  undisturbed  flow 
and  recovery  spectra  match,  indicating  that  the  effects  of  the  turbulence  pulse  are  purely 
transient.  Of  interest  are  the  very  large  spikes  in  the  spectra  of  the  recovery  period  plots 
for  y<  0.1 250  inches.  The  author  speculates  that  these  are  Tollmien-Schlichting  waves 
generated  during  the  laminar  to  turbulent  transition. 
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Figure  50.  Velocity  Spectra  and  Total  Power  Distribution,  50%  Chord,  y=  0.1250" 
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The  total  power  distribution  plots  show  trends  similar  to  the  velocity  spectra 
curves.  Near  the  surface  (Figure  46).  the  turbulent  power  is  quite  high  for  all  three  flow 
regimes.  Although  the  turbulent  pulse  does  not  affect  the  recovery  period  power  level, 
the  frequency  distribution  of  the  power  is  clearly  altered  by  the  turbulence  pulse.  The 
undisturbed  flow  curve  rises  slowly  to  the  asymptote  while  the  recovery'  period  curve 
rises  sharply  and  levels  off  very  quickly.  Near  the  middle  of  the  boundary  layer 
(Figure  48  and  Figure  49),  the  recovery  period  turbulence  power  is  much  lower  than 
that  of  the  free  stream.  This  shows  that  there  is  a  definite  recovery  process  operating  in 
this  region  as  the  flow  is  energized  back  to  the  free  stream  condition.  In  the  external- 
flow  region  (Figure  51),  the  undisturbed  flow  and  recovery  period  cumulative  power  are 
nearly  equal  indicating  that  the  effects  of  the  turbulence  pulse  are  purely  transitory. 

The  reverse-transition  to  laminar  flow  (variously  called  reversion,  inverse  or  re¬ 
verse  transition,  relaminarization,  and  laminarization)  observed  in  Figure  43  is  identical 
to  the  phenomenon  observed  by  Howard  [Ref.  6]  for  a  wing  immersed  in  propeller  wash. 
It  appears  that  the  turbulence  pulse  triggers  some  mechanism  in  the  boundary  layer  to 
laminarize  the  inherently  turbulent  boundary  layer  in  this  investigation  as  the  propeller 
wash  did  in  Howard's  work. 

Narasimha  and  Sreenivasan,  in  their  comprehensive  survey  of  relaminarizing 
flows,  identified  four  major  reversion  mechanisms:  reversion  by  dissipation;  reversion  in 
stably  stratified  flows;  reversion  in  highly  accelerated  flows;  and  relaminarization  by  ro¬ 
tation.  [Ref.  16] 

Laminarization  by  sudden  acceleration  of  the  flow  in  the  near  surface  region 
appears  to  be  the  mechanism  at  work  in  this  instance.  Analysis  of  the  turbulence  pulse 
regime  of  Figure  44  reveals  instantaneous  accelerations  of  8.000  to  10.000  ft.  sec2  are  felt 
by  the  flow  near  the  surface.  The  flow  appears  to  be  accelerated  both  by  the  turbulence 
pulse  itself  during  recovery  from  the  velocity  deficit  and  by  momentum  transfer  down 
through  the  boundary'  layer  (see  previous  section). 

One  of  the  most  widely  used  parameters  that  have  proven  to  predict  the  onset 
and  course  of  laminarization  of  a  turbulent  boundary  layer  is: 


v  dUe(x) 
L«!(x)  <<x 


(10) 


where  K  is  the  acceleration  parameter,  Ue(x)  is  the  edge  free  stream  velocity,  v  the 

kinematic  viscositv,  and  — - —  the  instantaneous  velocitv  derivative.  K,  as  used  in 
'  dx 
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Equation  10.  is  purely  a  free  stream  parameter  that  carries  no  information  at  all  about 
the  boundary  layer  itself.  Numerous  investigators  have  proposed  that  values  of  K 
greater  than  2  x  10-6  to  3.5  x  10~6  will  give  rise  to  laminarization  in  low  speed 
incompressible  flows.  [Ref.  16] 

In  this  investigation,  the  accelerated  flow  is  observed  in  the  boundary  layer  and 
in  the  external  free  stream  as  well.  For  this  reason,  the  author  chose  to  evaluate  K  lo¬ 
cally  within  the  boundary  layer  as  a  means  of  predicting  laminarization.  Applied  in  this 
manner,  K  becomes: 


K  = 


v  dL'(x.y) 
U2(x,y)  dx 


(11) 


where  L'(x,y)  is  the  local  ensemble  average  velocity  at  height  y  in  the  boundary  layer. 

Values  of  the  acceleration  parameter,  K,  range  from  2.2  x  IGF6  at  the  near¬ 
surface  region  to  1.7  x  1(E6  in  the  upper  regions  of  the  boundary  layer  at  50%  chord 
using  local  velocity  values.  The  value  calculated  near  the  surface  is  clearly  sufficient  to 
induce  laminarization  of  the  turbulent  boundary  layer  at  this  position.  It  is  interesting 
to  note  that  K  is  much  higher  at  30%  chord  where  it  remains  fairly  constant  at  4.4  x 
10~6  over  the  entire  boundary  layer  region;  but  K  does  not  play  a  significant  role  at  30°  o 
chord  other  than  hastening  the  return  to  naturally  laminar  flow. 

3.  70%  Chord 

The  ensemble  average  velocity  and  turbulence  profiles  plotted  in  Figure  52 
show  yet  entirely  different  characteristics  than  observed  at  30%  and  50To  chord.  The 
undisturbed  flow  is  fully  turbulent  with  a  high  level  of  turbulence  spanning  the  entire 
boundary  layer  region  (sample  200).  The  turbulence  pulse  flattens  the  velocity  profile 
while  inducing  high  turbulence  (6%  -  8%)  uniformly  over  the  vertical  region  (samples 
350  and  410).  During  recovery,  the  velocity  profiles  become  laminar  like  (sample  720) 
as  the  turbulence  decreases  in  the  middle  and  upper  regions  of  the  boundary  layer.  The 
turbulence  traces  appear  transitional  turbulent  in  the  recovery  regime.  Sample  850 
shows  the  flow  becoming  fully  turbulent  again  with  the  boundary  layer  height  increasing 
and  with  increasing  turbulence  in  the  lower  and  middle  section  of  the  boundary  layer. 
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Note  that  the  turbulence  level  near  the  surface  remains  high  and  fairly  constant 
throughout  the  period.  Clearly  the  laminarization  process  is  at  work  here  also,  but  the 
effects  are  greatly  reduced. 

The  smoothed  velocity  time  histories  plotted  in  Figure  53  show  that  the  velocity 
deficit  is  felt  throughout  the  boundary  layer.  The  velocity  deficit  grows  near  the  mid¬ 
layer  and  decreases  again  near  the  surface.  Of  interest  in  the  velocity  time  histories  is 
the  nature  of  the  velocity  increase  following  the  turbulence  pulse.  The  velocity  increase 
is  much  higher  at  70%  chord  than  at  the  previous  wing  locations.  Note  also  the  long 
recovery  period  and  slow  rise  and  decay  time  associated  with  the  velocity  increase. 

Figure  54  shows  the  smoothed  turbulence  time  histories.  Notice  the  high  tur¬ 
bulence  levels  that  persist  in  the  near-surface  regions  during  the  entire  period.  The  trace 
at  y  =  0.0702  inches  would  tend  to  indicate  that  the  turbulence  pulse  affects  the  mid¬ 
layer  region  more  than  the  near-surface  region  at  this  chord  location. 

The  velocity  spectra  and  total  power  distribution  graphs  are  shown  in 
Figure  55  through  Figure  60.  In  the  near-surface  region  (Figure  55  through 
Figure  57),  the  turbulence  pulse  mainly  affects  the  velocity  spectrum  of  the  recovery- 
period  while  total  power  remains  high  and  nearly  equal  over  the  three  flow  regimes. 
Toward  the  middle  of  the  layer  (Figure  58  and  Figure  59),  the  total  power  levels  of  the 
undisturbed  flow  and  recovery  period  regimes  drop  and  the  velocity  spectra  align  more 
closely.  It  is  interesting  to  note  that  the  total  power  of  the  recovery  period  is  less  than 
that  for  the  undisturbed  flow  in  the  middle  boundary  layer  region.  This  indicates  that 
there  is  a  definte  recovery  process  at  work  here.  In  the  external-flow  region 
(Figure  60).  the  turbulence  pulse  does  not  affect  the  velocity  spectrum  of  the  recovery- 
period.  Here  the  power  curves  for  the  undisturbed  How  and  recovery  period  regimes  are 
very  low  magnitude  and  equal. 

The  local  acceleration  parameter  (K)  from  Equation  11,  varies  from  1.3  x  10~6 
near  the  surface  to  1.9  x  10~6  at  the  middle  of  the  boundary  layer.  The  low  value  of  K 
near  the  surface  indicates  that  this  region  should  not  laminarize.  This  is  borne  out  by 
the  data;  the  near-surface  region  is  dominated  by  high  turbulence  and  the  mean  velocity 
responds  very  slowly  to  the  turbulence  pulse.  Conversely,  the  middle  boundary  layer 
region,  with  a  value  of  K  close  to  that  necessary  for  laminarization,  has  a  laminar-like 
response:  a  sharp  drop  in  turbulence  following  the  turbulence  pulse  with  a  corresponding 
sharp  velocity  increase. 

Apparently  the  turbulence  pulse  prompts  the  middle  boundary  layer  toward 
laminar  flow  but  the  pulse  is  unable  to  fully  penetrate  into  the  thick,  highly  turbulent 
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Sample  Number 


Figure  54.  Mean  Turbulence  Time  Histories  in  the  Boundary  Layer,  70 
Chord:  (Fs  =  50KHz) 
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Figure  57. 


Velocity  Spectra  and  Total  Pouer  Distribution,  70%  Chord,  y  =  0.0402" 
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Figure  59. 


Velocity  Spectra  and  Total  Poner  Distribution,  70%  Chord,  y  =  0.1941" 
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2000  4000  6000  8000  10000 

Frequancy  (Hz) 


near-surface  region.  Regardless,  the  turbulence  pulse  appears  able  to  prompt  a 
laminar-like  response  at  this  chord  location,  contrary  to  what  one  would  intuitively  ex¬ 
pect. 
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V.  CONCLUSIONS 


The  time-varying  effects  of  periodic  free  stream  turbulence  have  been  investigated 
on  a  wing  boundary  layer  at  three  representative  chord  locations  (laminar, 
transitional  turbulent,  and  turbulent)  on  the  wing.  An  experimental  method  was  devel¬ 
oped  and  tested  to  generate  non-thrusting,  uniform  periodic  turbulence  pulses  in  the 
external  flow  over  the  test  wing.  Although  the  turbulence  pulses  were  restricted  to  one 
test  condition,  50  pulses  per  second  and  a  turbulence  intensity  of  near  10°  o,  they  dem¬ 
onstrated  the  dynamic  character  of  the  boundary  layer  subjected  to  single-event  turbu¬ 
lence.  The  results  of  this  investigation  are  detailed  below. 

The  turbulence  pulse  was  a  short  (5  millisecond),  high  intensity  (10.5%  Tu),  low 
frequency  dominated  (no  significant  spectral  components  above  1  KHz),  non-thrusting 
event. 

The  turbulence  pulse  was  shown  to  fully  penetrate  the  boundary  layer  at  all  three 
chord  locations.  It  was  observed  that  the  turbulence  pulse  suppressed  the  high  fre¬ 
quency  spectral  components  in  the  recovery  period  regime.  Actual  penetration  was  less 
near  the  surface  in  the  turbulent  flow  region  of  the  wing,  but  the  effects  of  the  pulse  were 
still  quite  evident. 

The  boundary  layer  time-response  showed  a  distinct  transition  from  undisturbed 
(laminar,  transitional,  or  turbulent)  flow  to  highly  turbulent  flow  when  subjected  to  the 
turbulence  pulse  at  all  three  chord  locations. 

The  recovery  process  leading  back  to  undisturbed  flow  was  highly  dependent  on  the 

nature  of  the  undisturbed  flow  as  shown  below: 

Laminar:  Quick  boundary  layer  response  with  a  rapid  return  to  undisturbed  laminar 
flow. 

Transitional  Turbulent:  Slower  boundary  layer  response  than  in  the  laminar  region. 
Lanunarization  was  evident  during  the  recovery  process  followed  by  a  slow  return  to 
undisturbed  transitional  turbulent  flow.  The  boundary  layer  thickness  was  decreased 
during  the  recovery  period,  and  the  turbulence  intensity  dropped  to  1  6  of  its  normal 
undisturbed  value  near  the  surface. 

Turbulent:  Slow  boundary  layer  response.  Laminarization  is  quite  evident  in  the  mid 
and  upper  portions  of  the  boundary  layer  and  evident  to  a  lesser  degree  near  the  sur¬ 
face.  There  was  a  marked  decrease  in  boundary  layer  thickness  during  the  recovery 
period  in  this  region  before  returning  to  the  undisturbed  turbulent  value. 
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Laminarization  appears  to  be  mainly  prompted  by  rapid  acceleration  of  the  flow 
internal  to  the  boundary  layer.  The  internal  flow  is  accelerated  both  by  the  turbulence 
pulse  itself  during  recovery  from  the  velocity  deficit  and  by  momentum  transfer  down 
through  the  boundary  layer.  The  velocity  time-history  traces  show  that  the  small  ve¬ 
locity  deficit  associated  with  the  turbulence  pulse  in  external  flow  is  amplified  in  the 
middle  boundary  layer  region  due  to  momentum  transfer  down  through  the  layer.  The 
near-surface  flow  experiences  a  rapid  velocity  increase  as  a  result  of  the  momentum 
transfer. 

The  total  spectral  power  content  of  the  flow  regimes  (undisturbed,  flow,  turbulence 
pulse,  and  recovery  period)  was  more  dependent  on  vertical  position  in  the  boundary 
layer  than  on  chord  location.  Near  the  surface,  the  total  power  of  all  three  regimes  was 
very  high  and  nearly  equal  even  though  the  spectral  distribution  of  the  power  varied 
considerably.  The  middle  boundary  layer  region,  on  the  other  hand,  is  characterized  by 
a  high  power  turbulence  pulse,  medium  power  level  undisturbed  flow,  and  a  low  power 
recovery  period.  The  very  low  power  level  observed  in  the  recovery  period  shows  that 
the  turbulence  pulse  drives  the  flow  to  a  very  low  turbulent  energy  state.  Laminarization 
occurs  as  there  is  insufficient  spectral  power  to  support  turbulent  flow  at  this  point.  The 
flow  slowly  gains  power  and  transitions  back  to  the  undisturbed  flow  regime  as  the 
steady  state  characteristics  once  again  prevail.  At  the  edge  of  the  boundary  layer,  the 
power  level  of  the  turbulence  pulse  remains  quite  high  while  the  power  levels  of  the  un¬ 
disturbed  flow  and  recovery  regimes  are  low  and  nearly  equal. 
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VI.  RECOMMENDATIONS 


This  investigation  demonstrated  the  stabilizing  nature  of  periodic  turbulence  on  the 
transitional  and  turbulent  regions  of  the  wing.  The  study  was  done  at  a  single  low  angle 
of  attack  at  three  representative  chord  locations  (laminar,  transitional,  and  turbulent). 
Further  experimental  work  should  be  done  at  additional  chord  locations  and  at  higher 
angles  of  attack  to  further  study  the  stabilizing  effect  of  periodic  turbulence  on  attached 
boundary  layers  and  on  flow  near  separation. 

More  research  also  needs  to  be  undertaken  to  better  understand  the  mechanism  at 
work  in  the  laminarization  process.  Future  studies  should  investigate  the  individual 
pulse  passage  cycles  in  greater  detail  to  determine  if  ensemble-averaging  techniques  ad¬ 
equately  represent  the  highly  dynamic  laminarization  process,  or  if  statistical  methods 
or  a  single  pulse  passage  cycle  analysis  would  be  more  instructive. 

Data  collection  was  extremely  time  consuming  with  the  present  configuration.  The 
following  modifications  should  be  implemented  for  future  investigations: 

Store  and  process  all  data  in  binary  integer  format  rather  than  in  ASCII  real  format. 
This  will  reduce  the  mass  data  storage  requirement  by  a  factor  of  nine  and  also  greatly 
reduce  data  transfer  and  processing  time. 

Expand  the  random  access  memory  of  the  data  acquisition  computer  so  that  the  data 
collection  and  traverser  control  software  routines  can  be  managed  by  one  personal 
computer.  This  will  allow  for  full  automation  of  the  test  sequence. 

Use  a  hot-film  probe  instead  of  a  hot-wire  probe  in  the  NPS  low  speed  wind  tunnel 
to  avoid  contamination  induced  drift  problems. 

Decrease  the  digitizing  rate.  This  investigation  has  shown  that  little  useful  spectral 
information  is  contained  in  the  frequencies  above  10  KHz.  As  such,  the  sampling  rate 
could  easily  be  reduced  to  20  KHz  with  no  loss  of  meaningful  data.  This  step  alone 
will  reduce  mass  storage  by  60?  o  with  a  comparable  reduction  in  data  handling  and 
processing  time. 

Future  research  should  also  be  directed  toward  investigating  the  effect  of  periodic 
or  single-event  turbulence  on  an  airfoil  near  stall  or  in  post-stall  flight  with  regard  to  loss 
of  attached  flow  and  loss  of  control.  Questions  to  be  answered  include: 

What  is  the  effect  of  the  observed  stabilizing  phenomena  on  flow  near  separation? 

Can  the  results  be  worked  to  an  advantage  at  high  angles  of  attack  with  regard  to 
aircraft  supermancuverability  or  controllability?  This  is  important  as  augmented  lift 
from  dynamic-stall  maneuvering  will  have  limited  application  unless  the  pilot  can 
maintain  control  to  transition  quickly  from  one  maneuvering  state  to  another. 
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Directly  related  areas  of  interest  include  continued  studies  of  the  effects  of  wind 
gusts,  ship  airwake  turbulence,  missile  wake,  and  aircraft  wake  on  aircraft  stability  and 
controllability.  Future  studies  of  canard  vvake  airfoil  interaction  at  high  lift  seem  perti¬ 
nent.  especially  in  relation  to  the  continuing  supermaneuverability  studies  at  NFS. 
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APPENDIX  A.  AIRFOIL  COORDINATES 


UPPER 

x/c 

y/c 

x/c 

LOWER 

y/c 

1.  000 

.  0000 

1.  000 

.  0000 

.  990 

.  0008 

.  990 

0002 

.  980 

.  0015 

.  980 

0005 

.  970 

.  0022 

.  970 

0008 

.  960 

.  0030 

.  960 

-.  0010 

.  950 

.  0039 

.  950 

0011 

.  940 

.  0050 

.  940 

0011 

.  930 

.  0064 

.  930 

0009 

.  920 

.  0080 

.  920 

0007 

.  910 

.  0100 

.  910 

-. 0003 

.  900 

.  0121 

.  900 

.  0000 

.  890 

.  0142 

.  890 

.  0006 

.  850 

.  0220 

.  880 

.  0010 

.  800 

.  0310 

.  850 

.  0023 

.  750 

.  0400 

.  800 

.  0042 

.  700 

.  0483 

.  750 

.  0056 

.  650 

.  0561 

.  700 

.  0066 

.  600 

.  0635 

.  650 

.  0071 

.  550 

.  0704 

.  600 

.  0072 

.  500 

.  0766 

.  550 

.  0067 

.  450 

.  0816 

.  540 

.  0065 

.  400 

.  0858 

.  530 

.  0063 

.  350 

.  0894 

.  520 

.  0061 

.  300 

.  0918 

.  510 

.  0059 

.  250 

.  0927 

.  500 

.  0056 

.  240 

.  0926 

.  490 

.  0053 

.  230 

.  0923 

.  480 

.  0050 

.  220 

.  0919 

.  470 

.  0047 

.  210 

.  0914 

.460 

.  0044 

.  200 

.  0907 

.  450 

.  0040 

.  190 

.  0900 

.  440 

.  0036 

.  180 

.  0888 

.  430 

.  0031 

.  170 

.  0876 

.420 

.  0027 

.  160 

.0862 

.410 

.  0022 

.  155 

.  0854 

.  400 

.  0018 

.  150 

.0846 

.  390 

.  0013 

.  145 

.  0837 

.  380 

.  0009 

.  140 

.  0828 

.  370 

.  0004 

.  135 

.  0818 

.  360 

.  0000 

.  130 

.0808 

.  340 

0008 

.  125 

.  0797 

.  320 

0015 

.  120 

.  0785 

.  300 

-. 0020 

.  115 

.  0773 

.  280 

0025 

.  110 

.  0760 

.  260 

-. 0029 

.  105 

.  0746 

.  240 

-. 0032 

.  100 

.  0731 

.  220 

0034 
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095 

.  0716 

.  210 

0035 

090 

.  0700 

.  200 

0034 

085 

.  0684 

.  180 

0033 

080 

.  0666 

.  150 

0030 

075 

.  0647 

.  100 

0016 

070 

.  0628 

.  080 

0008 

060 

.  0585 

.  060 

0000 

050 

.  0537 

.  050 

.  0002 

040 

.  0481 

.  040 

.  0003 

030 

.  0415 

.  030 

.  0000 

020 

.  0335 

.  025 

0002 

015 

.  0286 

.  020 

0008 

010 

.  0233 

.  015 

0016 

005 

.  0169 

.  010 

0008 

.  005 

.  0020 
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APPENDIX  B.  DATA  SAMPLE  SPACING 


Sample  30%  Chord  50%  Chord  70%  Chord 

Number  (in)  (in)  (in) 


01 

0. 0032 

02 

0. 0033 

03 

0. 0038 

04 

0. 0046 

05 

0. 0057 

06 

0. 0071 

07 

0. 0087 

08 

0. 0107 

09 

0. 0130 

10 

0.  0156 

11 

0.  0184 

12 

0. 0217 

13 

0.  0251 

14 

0. 0288 

15 

0. 0330 

16 

0. 0372 

17 

0. 0418 

18 

0. 0467 

19 

0. 0518 

20 

0.  0572 

21 

0. 0628 

22 

0. 0687 

23 

0.  0748 

24 

0. 0811 

25 

0. 0877 

26 

0. 0944 

27 

0. 1014 

28 

0. 1087 

29 

0. 1162 

50 

0. 1238 

31 

0. 1316 

32 

0. 1396 

33 

0.  1477 

34 

0. 1560 

35 

0. 1644 

36 

0. 1731 

37 

0. 1818 

38 

0. 1906 

39 

0. 1996 

40 

0. 2087 

41 

42 

0. 0032 

0. 0048 

0. 0034 

0. 0060 

0. 0041 

0. 0075 

0.  0051 

0. 0095 

0. 0064 

0. 0122 

0. 0083 

0. 0154 

0. 0106 

0. 0192 

0. 0132 

0. 0236 

0.  0163 

0. 0286 

0. 0197 

0.  0342 

0. 0236 

0. 0427 

0. 0278 

G. 0469 

0. 0324 

0.  0541 

0. 0374 

0. 0619 

0. 0428 

0. 0702 

0. 0486 

0. 0790 

0. 0547 

0.  0884 

0. 0612 

0. 0982 

0. 0680 

0. 1086 

0.  0751 

0. 1194 

0.  0826 

0. 1307 

0. 0904 

0. 2425 

0. 0986 

0.  1547 

0. 1071 

0. 1674 

0. 1158 

0. 1805 

0. 1250 

0.  1941 

0. 1343 

0. 2080 

0. 1440 

0. 2223 

0. 1538 

0. 2371 

0. 1640 

0. 2521 

0. 1743 

0. 2676 

0.  1850 

0. 2833 

0. 1958 

0. 2994 

0. 2070 

0.  3157 

0. 2182 

0. 3324 

0. 2297 

0. 3493 

0. 2413 

0. 3665 

0. 2531 

0. 3839 

0.  2651 

0. 2722 

0. 2894 

0. 3017 

0.4016 
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APPENDIX  C.  PROGRAM  TRAVERSE 


100  REM 
110  REM 
120  REM 
130  REM 
140  REM 
150  REM 
160  REM 
170  REM 
180  REM 
190  REM 
200  REM 
210  REM 
220  REM 
230  REM 
240  REM 


*• 


•*** 


•*• 


PROGRAM  TRAVERSE.  BAS  ****************** 


THIS  PROGRAM,  WRITTEN  BY  D.  KINDELSPIRE  AND  MODIFIED  BY  R.  RENOUD 
(SEP  88)  CONTROLS  THE  VELMEX  8300  3-D.  TRAVERSER.  THE  INITIAL 
PORTION  OF  THE  PROGRAM  IS  DEDICATED  TO  OPERATOR-SELECTABLE  OPTIONS 
(I.E.  MOTOR  VELOCITY,  MOTOR  ACCELERATION,  ETC.)  AVAILABLE  FOR 
EACH  MOTOR.  DEFAULT  MOTOR  VALUES  FOR  EACH  OPTION  ARE  PROVIDED 
IN  THE  PROGRAMMING  AND  THE  OPERATOR  CAN  SELECT  THE  DEFAULTS  FOR 
ALL  MOTORS  WITH  ONE  KEYSTROKE.  THEW  PROGRAM  GIVES  THE  OPERATOR 
A  CHOICE  OF  MANUAL  INPUT  (MEANING  ONE  MOTOR  MOVEMENT  FOR  EACH 
INPUT  WITH  INPUT  IN  INCHES)  OR  COMPUTER -CONTROLLED  MOVEMENT 
THROUGH  THE  BOUNDARY  LAYER.  THE  SINE  FUNCTION-BASED  ALGORITHIM 
USED  TO  COMPUTE  THE  MOTOR  MOVEMENTS  THROUGH  THE  BOUNDARY  LAYER 
GIVES  MOVEMENT  COMMANDS  IN  INCREMENTS  OF  THE  MINIMUM  MOTOR 
MOVEMENT  STEP  SIZE  (0.000125  INCHES). 


250  REM 

260  RE  M  V“,“Vy'V**Wr7Wr‘V*,*«'w«  y- ■sW? 

1000  REM 

1010  REM  OPEN  THE  COM  PORT  AND  INITIALIZE  THE  MOTOR  SETTINGS 
1020  OPEN  "coral: 1200, n, 8,1, rs,cs,ds,cd"  FOR  RANDOM  AS  #1 
1030  REM  SET  MOTOR  DEFAULT  VALUES 

1040  DATA  1000,1000,1000,2,2,2,0.000125,0.000125,0.000125 
1050  READ  VI,  V2 ,  V3,  Rl,  R2,  R3,  Cl,  C2,  C3 

1060  DIM  Z(100),  Zl(lOO),  HWV(IOO),  RMS(IOO),  VEL(IOO),  NDVEL(IOO),  VRMS(IOO) 
1070  PI  =  3.  141592654?? 

1080  PRINT 

^  Q  ^  Q  P  R  X  \T  ^  ^  VwV  “V«V«V  "V“V  V»  »Wr*W*  #V“W*  *V  »V  V*  Vp*V«W»  »V  *  * 

USER  MUST  SELECT  'CAPS  LOCK'  FUNCTION  **" 

■-****************************************" 

DISPLAY  MOTOR  DEFAULT  SETTINGS 

INITIALIZED  VALUES  FOR  ALL  MOTOR  SETTINGS:" 

VELOCITY  =  1000  STEPS/SEC" 

RAMP( MOTOR  ACCELERATION)  =  2  (6000  STEPS/SEC 

DEFAULT  INCREMENTAL  UNITS  ARE  INCHES" 

* *****************  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 


1100  PRINT 
1110  PRINT 
1120  RLM 
1130  PRINT 
1140  PRINT 
1150  PRINT 
1160  PRINT 
1170  PRINT 


ii**,. 

"*v 

"*•; 

ii 

ii 

ii 


1180  PRINT  ***-. 

1190  PRINT 

1200  PRINT  "DO  YOU  WANT  TO  INPUT  YOUR  OWN  INCREMENT  EACH  TIME  OR  HAVE  THE" 


1210  PRINT 
1220  PRINT 
1230  PRINT 


COMPUTER  COMPUTE  MOVEMENTS  THROUGH  A  BOUNDARY  LAYER  FOR  YOU 


NOTE  *************" 


1240  PRINT  "FOR  COMPUTED  BOUNDARY  LAYER  MOVEMENT,  YOU  MUST  KNOW:" 


WHICH  MOTOR  WILL  MOVE  THE  APPROPRIATE  DIRECTION.  " 

HOW  THICK  (IN  INCHES)  IS  THE  B. L. " 

HOW  MANY  DATA  COLLECTION  POINTS  DO  YOU  WANT  IN  THE  B. L.  " 


1250  PRINT  1. 

1260  PRINT  "  2. 

1270  PRINT  "  3. 

1280  PRINT  " 

1290  PRINT 
1300  PRINT  "NOTE!! 

1310  PRINT  " 

1320  PRINT 

1330  INPUT  "MANUAL  CONTROL  OR  COMPUTER  CONTROL  (ENTER  ’MAN’  or  'CP')";  CON$ 


USE  MANUAL  CONTROL  TO  INITIALIZE  PROBE  POSITION  BEFORE' 
SELECTING  COMPUTER  CONTROLLED  MOVEMENT. 
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1340  IF  CON$  =  "CP"  THEN'  3130 

1350  REM  OPTION  TO  CHANGE  DEFAULT  SETTINGS  OF  VELOCITY  OR  ACCELERATION  RAMP 
1360  PRINT 
1370  PRINT 

13S0  PRINT  "  DO  YOU  WANT  TO  CHANGE  THE  VELOCITY  OR  ACCELERATION  RAMP" 

1390  PRINT  "  DEFAULT  SETTINGS?  (Y  or  N)" 

1400  PRINT 

1410  PRINT  "IF  ’NO',  THIS  PROGRAM  WILL  THEN  LET  YOU  DEFINE  THE" 

1420  PRINT  "DISTANCE  YOU  WANT  TO  MOVE  (IN  INCHES).  IF  'YES'," 

1430  PRINT  "YOU  CAN  CHANGE  ANY  OR  ALL  OF  THE  DEFAULT  SETTINGS  FOR  ANY  MOTOR.  " 
1440  PRINT 
1450  PRINT 
1460  PRINT 

1470  INPUT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THE  DEFAULT  SETTINGS?  (Y  or  N)";  D$ 
1480  IF  D$  =  "Y"  THEN  1520 
1490  IF  D$  =  "N"  THEN  2150 
1500  REM 

1510  REM  ****  OPERATOR  SELECTED  MOTOR  VARIABLES  ***** 

1520  PRINT 
1530  PRINT 

1540  INPUT  "WHICH  DEFAULT  VALUE?  (ENTER  ' 1 ’ FOR  VELOC  OR  '2'  FOR  ACCEL  RAMP)";  L 
1550  ON  L  GOTO  1620,  1860 

1560  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  VELOCITY’?  (Y  OR  N)" 

1570  INPUT  V$ 

1580  IF  VS  =  "Y"  THEN  1620 

1590  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  ACCELERATION  RAMP?  (Y  or  N)" 

1600  IF  RS  =  "Y"  THEN  1920 
1610  IF  RS  =  "N"  THEN  1380 
1620  PRINT 
1630  PRINT 

1640  INPUT  "WHICH  MOTOR  VELOCITY  DO  YOU  WISH  TO  CHANGE?  (1,2,  or  3)";  J 
1650  ON  J  GOTO  1660,  1760,  1810 
1660  PRINT 
1670  PRINT 

168U  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  #l";  VI 
1690  PRINT 
1700  PRINT 
1710  PRINT 

1720  PRINT  "DO  YOU  WANT  TO  CHANGE  VELOCITY  OF  ANOTHER  MOTOR?  (Y  OR  N)" 

1730  INPUT  VS 

1740  IF  VS  =  "Y"  THEN  1620 
1750  IF  VS  =  "N"  THEN  1360 
1760  PRINT 
1770  PRINT 

1780  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  2";  V2 

1790  PRINT 

1800  GOTO  1710 

1810  PRINT 

1820  PRINT 

1830  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  #3";  V3 

1840  PRINT 

1850  GOTO  1710 

1860  PRINT 

1870  PRINT 

1880  INPUT  "WHICH  MOTOR  ACCEL  RAMP  DO  YOU  WANT  TO  CHANGE?  (1,  2,  or  3)";  K 
1890  ON  K  GOTO  1900,  1990,  2050 
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1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 


PRINT 

PRINT 

INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #l";  R1 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  THE  ACCEL  RAMP  OF  ANOTHER  MOTOR?  (Y  or  N)?" 
INPUT  RM$ 

IF  RMS  =  "Y"  THEN  1860 
IF  RM$  =  "N"  THEN  1380 
PRINT 
PRINT 

INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #2";  R2 

PRINT 

PRINT 

GOTO  1930 

PRINT 

PRINT 

INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #3";  R3 

PRINT 

PRINT 

GOTO  1930 

REM 

REM  DEFINE  DISTANCE  TO  MOVE  MOTOR 

PRINT 

PRINT 

PRINT 

REM  INITIALIZE  MOTOR  INCREMENTS  TO  ZERO 

11  =  0 

12  =  0 
13  a  0 
PRINT 

PRINT  "  *************************************************************' 


PRINT  "  ** 
PRINT  "  ** 
PRINT  "  ** 
PRINT  "  ** 
PRINT  "  ** 
PRINT  "  ** 


DEFINE  WHICH  MOTOR  YOU  WANT  TO  MOVE 

NOTE!!!  A  POSITIVE  (’+')  INCREMENT  TO  A  MOTOR 
MOVES  TRAVERSER  AWAY  FROM  THAT  PARTICULAR  MOTOR 


PRINT  "  ** 
PRINT  "  ** 


--  MOTOR  #1  MOVES  THE  PROBE  UPSTREAM  AGAINST  THE  FLOW  ** 
--  MOTOR  in  MOVES  THE  PROBE  TOWARD  THE  ACCESS  WINDOW 
--  MOTOR  MOVES  THE  PROBE  VERTICALLY  DOWNWARD 
PRINT  "  ************************************************************* 
PRINT 
PRINT 

INPUT  "WHICH  MOTOR  DO  YOU  WANT  TO  MOVE?  (1,2,  or  3)";  L 

ON  L  GOTO  2350,  2460,  2530 

PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #1?" 

PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

INPUT  II 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

INPUT  C$ 

IF  C$  =  "Y"  THEN  2300 
IF  C$  =  "N"  THEN  2600 
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2460  PRINT 
2470  PRINT 

2480  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #2?" 

2490  PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

2500  INPUT  12 
2510  PRINT 
2520  GOTO  2410 
2530  PRINT 
2540  PRINT 

2550  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #3?" 

2560  PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

2570  INPUT  13 
2580  PRINT 
2590  GOTO  2400 
2600  PRINT 
2610  PRINT 

2620  REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 
2630  PRINT  "  *********************************" 

2640  PRINT 

2650  PRINT  "SUMMARY  OF  OPERATOR  INPUTS:  " 

2660  PRINT  "  MOTOR  #1  VELOCITY  =  VI 

2670  PRINT  "  ACCELERATION  RAMP  =  R1 

2680  PRINT  "  INCREMENTAL  DISTANCE  =  II;  "INCHES" 

2690  PRINT  "  MOTOR  #2  VELOCITY  =  V2 

2700  PRINT  "  ACCELERATION  RAMP  =  R2 

2710  PRINT  "  INCREMENTAL  DISTANCE  =  12;  "INCHES" 

2720  PRINT  "  MOTOR  #3  VELOCITY  =  V3 

2730  PRINT  "  ACCELERATION  RAMP  =  R3 

2740  PRINT  "  INCREMENTAL  DISTANCE  =  13;  "INCHES" 

2750  PRINT 

2760  PRINT  "  ********************************* 1  f 
2770  PRINT 
2780  PRINT 

2790  PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N)" 

2800  PRINT 

2810  PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.  ENTER  'y'  TO  RETURN" 

2820  PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

2830  INPUT  V$ 

2840  IF  V$  =  "Y"  THEN  1360 
2850  GOSUB  3060 
2860  PRINT 
2870  PRINT 

2880  INPUT  "DO  YOU  WANT  TO  INPUT  ANOTHER  MANUAL  MOTOR  MOVEMENT  (Y  or  N)";  M$ 
2890  IF  M$  =  "Y"  THEN  2140 
2900  PRINT 

2910  PRINT  "DO  YOU  WANT  TO  INPUT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT?" 

2920  PRINT  "  *********  NOTE! !  !  *********  " 

2930  PRINT  "  ALL  PREVIOUS  MOTOR  INCREMENT  INPUTS  HAVE  BEEN  ZEROIZED." 

2940  PRINT  "PROGAM  WILL  LET  YOU  CHOOSE  MANUAL  OR  CP-CONTROLLED  B.  L.  MOVEMENT. 
2950  PRINT  "*****  (IF  'NO',  THE  PROGRAM  WILL  END).  *****" 

2960  PRINT 

2970  INPUT  "DO  YOU  WANT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT  (Y  or  N)";  N$ 

2980  IF  N$  =  "Y"  THEN  1220 
2990  PRINT 
3000  PRINT 
3010  PRINT 
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PR  I  NT  ' '  **************************** 1 ' 

PRINT  "  THE  PROGRAM  HAS  ENDED. " 

PR  I  NT  "  ****************************" 

END 

REM  *******  MOTOR  MOVEMENT  SUBROUTINE  ********* 

PRINT  #1,  "fit":  PRINT  #1 ,  "E";  "Cl=";  Cl;  ":C2=";  C2;  C3=";  C3 

PRINT  #1,  ”I1=";  II;  ":V1=";  VI;  ":R1=":  Rl; 

PRINT  #1,  ":I2=";  12;  ":V2=";  V2;  ":R2=;  R2 
PRINT  #1,  "I3=";  13;  ":V3=";  V3;  ":R3=";  R3; 

RETURN 

REM  ******************************************** 

PRINT 

REM  *******  COMPUTER  CONTROLLED  BOUNDARY  LAYER  MOVEMENT  ******* 

PRINT 

PRINT 

PRINT  "  ****************************************" 

PRINT  "  **  NOTE  !  !  !  **" 

PRINT  "  **  COMPUTER  CONTROLLED  BOUNDARY  LAYER  **" 

PRINT  "  **  MOVEMENT  IS  PROGRAMMED  WITH  A  **" 

PRINT  "  **  DEFAULTED  NEGATIVE  MOTOR  INCREMENT  **" 

PRINT"  **  (i. e.  MOTOR  #3  WILL  MOVE  UPWARD 

PRINT  "  **  BY  ENTERING  A  (+)  B. L.  DISTANCE).  **" 

PRINT  "  ************ **************************** 1 ' 

PRINT 

REM  SET  INITIAL  BL  HEIGHT  AND  NUMBER  OF  DATA  POINTS  TO  ZERO 
N  =  0 
BL  =  0 
PRINT 

INPUT  "ENTER  ESTIMATED  INITIAL  POINT  DISTANCE  ABOVE  SURFACE  (INCHES)";  ZIP 
PRINT 

INPUT  "WHAT  IS  THE  B.  L.  THICKNESS  (IN  INCHES)  THAT  YOU  WANT  TO  MEASURE";  BL 
PRINT 

INPUT  "HOW  MANY  DATA  POINTS  DO  YOU  WANT  TO  SAMPLE  IN  THE  B.L.  ";  N 
PRINT 

INPUT  "WHICH  MOTOR  WILL  MOVE  YOU  THROUGH  THE  B.L.  (1,2,  or  3)";  MN 

ON  MN  GOTO  3380,  3730,  4090 

PRINT 

PRINT 

REM  MOTOR  1  CP -CONTROLLED  MOTOR  MOVEMENT 
Z(0)  =  0 

11  =  0 

12  =  0 
13  =  0 

FOR  J  =  1  TO  N 

Z(J)  =  BL  *  (1  -  SIN( ( PI  /  2)  *  (1  +  (J  -  1)  /  (N  -  1)))) 

STP  =  CINT( ( Z( J)  -  Z(J  -  1))  /  .000125)  *  .000125 

Z(J)  =  Z(J  -  1)  +  STP 

II  =  -STP 

PRINT 

PRINT 

PRINT  "  ************************" 

PRINT 

PRINT  "  DATA  POINT  ";  J 

PRINT 

GOSUB  3060 

PRINT 
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3580  PRINT 

3590  IF  J  =  N  THEN  3620 

3600  INPUT  "PRESS  'ENTER'  FOR  NEXT  MOVEMENT";  MOVE$ 

3610  NEXT  J 

3620  PRINT  "ALL  MOVEMENTS  COMPLETE" 

3630  PRINT 
3640  PRINT 

3650  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?" 

3660  PRINT 

3670  PRINT  "IF  'Y' ,  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

3680  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

3690  PRINT 

3700  INPUT  "ANOTHER  MOVEMENT";  R$ 

3710  IF  R$  =  "Y"  THEN  1290 

3720  IF  R$  =  "N"  THEN  2990 

3730  PRINT 
3740  PRINT 

3750  REM  CP-CONTROLLED  MOTOR  MOVEMENT  FOR  MOTOR  2 

3760  Z(0)  =  0 

3770  II  =  0 

3780  12  =  0 

3790  13  =  0 

3800  FOR  J  =  1  TO  N 

3810  Z(J)  =  BL  *  (1  -  SIN( (PI  /  2)  *  (1  +  (J  -  1)  /  (N  -  1)))) 

3820  STP  =  CINT( ( Z ( J)  -  Z(J  -  1))  /  .000125)  *  .000125 

3830  Z(J)  =  Z(J  -  1)  +  STP 

3840  12  =  -STP 

3850  PRINT 

3860  PRINT 

3870  PRINT 

3880  PRINT  "  ************************ ' ' 

3890  PRINT 

3900  PRINT  "  DATA  POINT  ";  J 

3910  PRINT 

3920  GOSUB  3060 

3930  PRINT 

3940  PRINT 

3950  IF  J  =  N  THEN  3980 

3960  INPUT  "PRESS  ENTER  FOR  NEXT  MOVEMENT";  MOVE$ 

3970  NEXT  J 

3980  PRINT  "ALL  MOVEMENTS  COMPLETED" 

3990  PRINT 
4000  PRINT 

4010  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?";  R$ 
4020  PRINT 

4030  PRINT  "IF  'Y' ,  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

4040  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

4050  PRINT 

4060  INPUT  "ANOTHER  MOVEMENT";  R$ 

4070  IF  R$  =  "Y"  THEN  1290 
4080  IF  R$  =  "N"  THEN  2990 
4090  PRINT 
4100  PRINT 

4110  REM  CP -CONTROLLED  MOTOR  MOVEMENT  FOR  MOTOR  3 
4120  Z(0)  =  0 
4130  II  =  0 
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4140  12  =  0 
4150  13  =  0 
4160  FOR  J  =  1  TO  N 

4170  Z(J)  =  BL  *  (1  -  SIN( (PI  /  2)  *  (1  +  (J  -  1)  /  (N  -  1)))) 

4180  STP  =  CINT((Z( J)  -  Z(J  -  D)  /  .000125)  *  .000125 

4190  Z(J)  =  Z(J  -  1)  +  STP 

4200  13  =  -STP 

4210  PRINT 

4220  PRINT 

4230  PRINT 

4240  PRINT  "  **********★■**********■***" 

4250  PRINT 

4260  PRINT  "  DATA  POINT  J 

4270  PRINT 

4280  GOSUB  3060 

4290  PRINT 

4300  PRINT 

4310  IF  J  =  N  THEN  4340 

4320  INPUT  "PRESS  'ENTER'  FOR  NEXT  MOTOR  MOVEMENT";  MOVE$ 

4330  NEXT  J 

4340  PRINT  "ALL  MOVEMENTS  COMPLETED" 

4350  PRINT 
4360  PRINT 

4370  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N) ?  ;  R$ 
4380  PRINT 

4390  PRINT  "IF  'Y',  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

4400  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

4410  PRINT 

4420  INPUT  "ANOTHER  MOVEMENT";  R$ 

4430  IF  R$  =  "Y"  THEN  1290 
4440  IF  R$  =  "N"  THEN  2990 
4450  END 


102 


APPENDIX  D.  PROGRAM  HWCAL 


100  REM 
110  REM 
120  REM 
130  REM 
140  REM 
150  REM 
160  REM 
170  REM 
180  REM 
190  REM 
200  REM 
210  REM 
220  REM 
230  REM 


yfiwV'VV'Vw'VfVfyr  VwV  PRO  GRAM  HW  CAL  S  A  S  jV*v^fTV*'*TWwf 

This  program  is  used  to  calibrate  the  hot  wire  in  the  low  speed  wind 
tunnel.  The  hotwire  is  calibrated  against  a  pitot/static  probe 
located  in  the  free  stream  adjacent  to  the  hotwire.  Pitot/static 
output  (delta  P)  is  in  cm  of  water.  Hotwire  output  from  the  IFA-100 
is  digitally  read  by  a  Metrabyte  16f  A  to  D  board  controlled  by 
Labtech  Notebook  software.  Suggested  A  to  D  speed  is  1000  Hz  with  a 
period  of  1.  01  sec.  The  first  1000  data  points  are  used  by  this 
program.  The  program  output  file  is  directly  compatable  with 
"Grapher"  for  plotting  calibration  curves  and  determining 
conversion  equations. 


250  REM 

260  REM  VARIABLES: 

270  REM  dp  :  Pitot/static  probe  delta  P  (cm  water) 

280  REM  hwv  :  Hotwire  voltage  (volts) 

290  REM  vfps:  Tunnel  velocity  (feet  per  second) 

300  REM  ofi  :  Output  file  name 

310  REM 

320  REM  o^fV'V'jV^yfiVVfV'VrVf'VVwfVwVfwf^ffwV'VVfyfVfVfTW'Vf'Wc^V'yf^rf^'VfVf rfAVfVvVf^HfVf'WwWfjV'WoV'VV'V'y'VfiV 

1000  DIM  DP( 20) ,  HWV( 20) 

1010  R  =  53. 3:  REM  gas  constant  for  air 
1020  G  =  32. 174 

1030  CLS  :  PRINT  "Enter  output  file  name":  INPUT  ;  0FI$ 

1040  OPEN  0FI$  FOR  OUTPUT  AS  # 2 

1050  PRINT  n,  "’em  water"’,  "'hwv'",  "' vi%km"  .5’" 

1060  1=1 

1070  PRINT  :  PRINT  "Enter  cm  of  water  for  calibration  point  I 
1080  INPUT  ;  DP(I) 

1090  PRINT  :  HWV(I)  =0 

1100  REM  Execute  A  to  D  conversion  of  hotwire  using  Labtech  Notebook  "go" 
1110  SHELL  "  nb  go" 

1120  REM  Labtech  Notebook  writes  output  file  to  data. prn 

1130  OPEN  "c:  nb  data  data. prn"  FOR  INPUT  AS  #1 

1140  FOR  J  =  1  TO  1000 

1150  INPUT  #1,  HWV 

1160  HWV(I)  =  HWV(I)  +  HWV 

1170  NEXT  J 

1180  CLOSE  in 

1190  HWV(I)  =  HWV(I)  /  1000 

1200  CLS  :  PRINT  "hvw  for  point  I;  "  =  HWV(I) 

1210  INPUT  "another  calibration  point  (y/n)";  C$ 

1220  IF  C$  =  "n"  THEN  GOTO  1250 
1230  1=1+1 
1240  GOTO  1070 

1250  CLS  :  INPUT  "Enter  ambient  temperature  (Deg  F)";  T 
1260  T  =  T  +  460 
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1270  INPUT  "Enter  barometric  pressure  (in.  Hg)";  P 
1280  PRINT  "Wait  -  Writing  Output  Data  File" 

1290  P  =  P  /  29.92  *  14.7  *  144 
1300  R0  =  P  /  (R  *  T  *  G) 

1310  FOR  J  =  1  TO  I 

1320  VFPS  =  SQR( 2  *  DP(J)  *  2.048  /  RO) 

1330  PRINT  #2,  DP( J) ,  HWV(J),  VFPS,  HWV(J)  2,  SQR(VFPS) 
1340  NEXT  J 
1350  CLOSE  #2 

1360  CLS  :  INPUT  "Another  calibration  run  (y/n)  ;  C$ 
1370  IF  C$  =  "y"  THEN  GOTO  1000 
1380  PRINT  "BYE" 

1390  END 
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APPENDIX  E.  PROGRAM  BLDATA 


100 

110 

111 

112 

113 

140 

150 

160 

170 

180 

190 

200 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1125 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 


REM  ******* 
REM 


PROGRAM  BLDATA.  BAS  ******************** 


■VrVfVrVrVi WrtWWirsWIr 


REM  This  driver  program  controls  the  analog-to-digital  data  acquisition 
REM  process  using  the  software  package  LABTECH  NOTEBOOK  in  conjunction 
REM  with  the  Metrabyte  16f  A  to  D  card. 

REM  This  progar m  allowes  the  user  to  enter  the  number  of  data  points 
REM  for  the  run.  It  initiates  the  data  collection  by  running  LABTECH 

REM  NOTEBOOK  program  "go"  directly  from  the  nb  directory  via  dos 

REM  command.  This  control  program  also  generates  incremental  names  for 
REM  the  data  files. 

REM 

CLS 

INPUT  "Enter  number  of  data  points  for  this  wing  position  (1-99)";  N 
IF  (N  <  1)  OR  (N  >  99)  GOTO  1010 

REM  ***  generating  string  segments  for  data  file  names 
BS  =  MID$( STR$( 1 ) ,  2):  REM  string  number  "l" 

E$  =  MIDS( STR$(N) ,  2):  REM  ending  string  number  "n" 

X$  =  "XXXXXX" 

EX$  =  ".DAT" 

CLS 

PRINT  "Data  files  will  be  incremented  from:" 

PRINT  (X$  +  B$  +  EX$);  "  To  (X$  +  E$  +  EX$) 

PRINT 

PRINT  "Enter  common  data  file  id  (XXXXXX)  for  this  run 
PRINT  "(Six  characters  maximum  -  no  extension" 

INPUT  FS 

IF  LEN(F$)  >  6  OR  LEN(F$)  <  1  GOTO  1120 
CLS 

PRINT  N;  "Data  files  will  be  generated  and  incremented  as  follows:" 
PRINT  (F$  +  B$  +  EX$);  "  To  ";  (F$  +  E$  +  EX$) 

PRINT 

INPUT  "Are  the  number  of  points  and  file  names  ok  (y  n)";  C$ 

IF  CS  =  "n"  GOTO  1000 
IF  C$  =  "y"  GOTO  1230 
GOTO  1190 
CLS 

FOR  I  =  1  TO  N 

PRINT  :  PRINT  "ready  for  point  I;  "  of  "j  N;  "  points":  PRINT 
PRINT  "Insert  New  Data  Disk  in  Drive  A" 

PRINT  "Strike  any  key  when  ready" 

BEEP 

C$  =  INKEY$:  IF  C$  =  ""  GOTO  1290 

REM  ***  generating  incremented  data  file  name 

IF  (I  >  10)  OR  (I  =  10)  THEN  1$  =  MID$(STR$( I) ,  2) 

IF  (I  <  10)  THEN  1$  =  (MID$(STR$(0) ,  2)  +  MID$(STR$( I) ,  2)) 

FI$  =  (F$  +  1$  +  EX$) 

PRINT  "Doing  A  to  D  conversion  and  writing  data  file  ";  FI$ 

SHELL  "cd  nb":  REM  DOS  command  to  ensure  system  in  nb  subdirectory 
SHELL  "go":  REM  executing  labtech  notebook  via  DOS  from  NB  subdir 
Rl$  =  "ren  A:  DATA.  PRN  " 
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1380  R2$  =  Rl$  +  FI$ 

1390  SHELL  R2$ 

1400  NEXT  I 

1410  INPUT  "Another  run"; 
1420  IF  C$  =  "y"  GOTO  1000 
1430  PRINT  "BYE" 

1440  END 


APPENDIX  F.  PROGRAM  CONVERT 

)  PROGRAM  CONVERT 

C ********************************************************************** 

c 

C  THIS  PROGRAM  CONVERTS  THE  RAW  HOT-WIRE  VOLTAGE  AND  POSITION  SENSOR 
’  C  VOLTAGE  INTO  VELOCITY  DATA  AN  A  SYNCHRONIZATION  PULSE.  INPUT  DATA 

C  IS  IN  ASCII  INTEGER  FORM  WITH  THE  INTEGER  VALUES  0  TO  4095  REP- 
C  RESENTING  0  TO  5  VOLTS.  THE  BAR  SENSOR  VOLTAGE  WAS  PROCESSED  TO 
C  RETURN  A  VALUE  OF  2  AT  THE  EXACT  CENTER  OF  THE  POSTION  SENSOR 
C  PULSE  AND  A  0  ELSEWHERE.  THE  HOT-WIRE  VOLTAGE  TO  VELOCITY  CONVER- 
C  SION  CONSTANTS  ARE  OBTAINED  FROM  THE  HOT-WIRE  CALIBRATION  DATA 
C  FROM  PROGRAM  HWCAL. BAS. 

C 

C 

REAL  HWV,  PROP,  A,  B,  PPROP,  TIC 
CHARACTER* 12  INFILE,  OUTFILE,  CHOICE*2 
INTEGER  I , J,K 
DIMENSION  DAT  (50000,2) 

A  =  0 
B  =  0 
C 

5  PRINT  '(/A  )','  ENTER  THE  INPUT  DATA  FILE  NAME:  ' 

READ  (*, ' ( A) ' )  INFILE 
C 

PRINT  '(/A  )',  '  ENTER  THE  OUTPUT  DATA  FILE  NAME: ’ 

READ  ( * , ' ( A ) ' )  OUTFILE 
C 

1  OPEN  (10,  FILE=INFILE,  STATUS  =  'OLD') 

OPEN  (11,  FILE=OUTFILE,  STATUS  =  'NEW') 

C 

PRINT  '(/A  )' ,  '  CONVERSION  FROM  HWV  TO  VELOCITY  IS  OF  THE  FORM: ' 
PRINT  '(/A  )' ,  '  VEL  =  (A  *  HWV**2  +  B)  **  2' 

PRINT  '(/A/)' ,  '  .  .  .  ' 

WRITE  (*,*)  '  A  =  ' ,  A 
WRITE  (*,*)  '  B  =  ' ,  B 
PRINT  '(/A/)',  ’  CHANGE  A  OR  B  (Y/N)  ?  ’ 

READ  (*, 1 (A) ’ )  CHOICE 
IF  (CHOICE  .EQ.  'N')  GOTO  7 
PRINT  '(/A  )’ ,  '  A  =  ' 

READ  (*,*)  A 

PRINT  ’(/A  )' ,  '  B  =  ' 

READ  (*,*)  B 
C 

7  PPROP=0 

PRINT  '(/A/)’,  '  READING  DATA  INTO  ARRAY  -  WAIT' 

.  DO  10  J  =  1,  50000 

READ  (10,  *,  END=50 ,  ERR=50)  HWV,  PROP 

DAT  (J,l)  =  HWV 
DAT  ( J, 2)  =  PROP 
10  CONTINUE 
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c 

50  CONTINUE 

PRINT  '(/A/)',  '  DOING  DATA  CONVERSION  -  WAIT  SOME  MORE' 

DO  20  I  =  1,  J-l 

DAT  (1,1)  =  (A  *  (DAT  (1,1)  *  5.  /  4095.)  **  2  +  B)  **  2 
PROP  =  DAT  (1,2)  *  5.  /  4095. 

TIC  =  0. 

IF  (PROP  .  GT.  1.5  .AND.  PROP  .  LT.  2.  .AND.  PROP  .  LT.  PPROP) 
+TIC  =  2. 

PPROP  =  PROP 
DAT  (1,2)  =  TIC 
20  CONTINUE 

C 

PRINT  '(/A/)’,  '  WRITING  OUTPUT  FILE  -  ALMOST  DONE! 

DO  30  I  =  1,  J-l 

WRITE  (11,80)  DAT  (1,1),  DAT  (1,2) 

80  FORMAT  (F11.4,  F5.  0) 

30  CONTINUE 

C 

WRITE  (*,*)  '  #  OF  CONVERSIONS  IS  :  ,  (J-l)  *  2 
CLOSE  (10) 

CLOSE  (11) 

WRITE  (*.*)  '  ANOTHER  RUN  (Y/N)  ?' 

READ  (*,  ( A) ' )  CHOICE 
IF  (CHOICE  .EQ.  * Y* )  GOTO  5 
C 

END 


< 


108 


LIST  OF  REFERENCES 


1.  Mueller,  T.J.,  Low  Reynolds  S  umber  Vehicles,  AGARDograph  2SS, 

^  AGARD-AG-2S8,  pp.  1-16,  February-  19S5. 

2.  Meier,  H.U.,  and  Kreplin,  H.P.,  "Influence  of  Freestream  Turbulence  on 
Boundary- Layer  Development,"  AlAA  Journal,  Vol.  18,  No.  7,  pp.  11-15,  January 
19S0. 

3.  National  Aeronautics  and  Space  Administration  Final  Report  (NASA-CR-18063S). 
The  In  fluence  of  Tree-Stream  Turbulence  on  Separation  of  Turbulent  Boundary  Lay¬ 
ers  in  Incompressible,  Two-Dimensional  Flow,  by  J.L.  Potter,  and  others,  pp.  1-42. 
22  December  19S6. 

4.  Kindelspirc.  DAW.  The  Effects  of  Freestream  Turbulence  on  Airfoil  Boundary  Layer 

*  Behavior  at  Low  Reynolds  Xumbcrs,  Master's  Thesis,  Naval  Postgraduate  School. 

Monterey.  CA,  September  19S8. 

5.  Brendel.  M..  Experimental  Study  of  the  Boundary  Layer  on  a  Low  Reynolds  X  umber 
Airfoil  in  Steady  and  Unsteady  Flow,  Ph.D.  Dissertation.  University  of  Notre  Dame, 
Notre  Dame.  IN.  May  1986. 

6.  Floward.  Richard  M.,  An  Investigation  of  the  Effects  of  the  Propeller  Slipstream  on 
a  Wing  Boundary  Layer,  Ph.D.  Dissertation,  Texas  A&M  University,  College  Sta¬ 
tion,  TX,  May  1987. 

7.  Department  of  Aeronautics  and  Astronautics,  Laboratory  Manual  for  Low-Speed 
Wind-Tunnel  Testing,  Naval  Postgraduate  School,  Monterey,  CA,  October  1983. 

9 

8.  TSI  Incorporated,  IFA  100  System  Instruction  Manual,  Revision  C.  St.  Paul,  MN, 
August  1987. 


109 


9.  Brendat.  J.S.  and  Piersol,  A.G..  Engineering  Applications  of  Correlation  and  Spectral 
Analysis,  John  Wiley  and  Sons.  New  York..  NY.  19S0. 

10.  Johnson,  D.K.,  A  Data  Analysis  System  for  Unsteady  Turbulence  Measurements, 
Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  September  19S8. 

11.  Vast,  J.D.,  and  Fernholz,  H.H.,  "A  Discussion  of  Probe  Effects  and  Improved 
Measuring  Techniques  in  the  Near-Wall  Region  of  an  Imcompressible  Three- 
Dimensional  Turbulent  Boundary  Layer,"  Turbulent  Boundary  Layers  -  Exper¬ 
iments,  Theory  and  Modelling,  AGARD  Conference  Proceedings  No.  271,  pp.  10-1 
to  10-17,  January,  1980. 

12.  Roane.  D.P..  Jr.,  The  Effect  of  a  Turbulent  Airstream  on  a  Vertically  Launched 
Missile  at  High  Angles  of  Attack,  Master's  Thesis,  Naval  Postgraduate  School, 
Monterey,  CA.  December  19S7. 

13.  Thermo-Systems  Inc..  Hot  Hire  -  Hot  Film  -  Ion  Anemometer  Systems. 
CAT  FORM  6560375.  St.  Paul.  MN.  1975. 

14.  American  Institute  of  Aeronautics  and  Astronautics,  AIAA-S7-1259.  Influence  of 
a  Heated  Leading  Edge  on  Boundary  Layer  Growth,  Stability,  and  Transition,  by 
Landrum.  D.B.,  and  Macha.  J.M.,  pp.  1-9.  June  19S7. 

15.  Bradshaw,  P.,  An  Introduction  to  Turbulence  and  its  Measurement,  Pergamon  Press, 
New  York.  NY.  1971. 

16.  Narasimha,  R.,  and  Sreenivasan.  K.R.,  "Relaminarization  of  Fluid  Flows,"  Ad¬ 
vances  in  Applied  Mechanics.  Vol.  39,  pp.  246-256,  Academic  Press,  Inc.,  New  York, 
NY,  1979. 


110 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1.  Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  YA  22304-6145 

2.  Library.  Code  0142  2 

Naval  Postgraduate  School 

Monterey.  CA  93943-5002 

3.  Chairman  2 

Department  of  Aeronautics  and  Astronautics,  Code  67 

Naval  Postgraduate  School 
Monterey.  CA  93943-5000 

4.  Commander  1 

Naval  Air  Systems  Command 

Washington.  DC  20360 

5.  Commandant  (G-PTE-1 )  2 

United  States  Coast  Guard 

Washington.  DC  20593-0001 

6.  NASA  Langley  Research  Center  1 

MS  1S5  Technical  Library 

Hampton.  VA  23665 

7.  NASA  Ames  Research  Center  1 

Technical  Library 

Moffett  Field.  CA  94035 

S.  Dr.  Bruce  J.  Holmes  1 

Head.  Flight  Applications  Branch 
NASA  Lanclev  Research  Center 
MS  247 

Hampton,  VA  23665 

9.  Dr.  S.  J.  Milev  1 

c  o  Aerospace  Engineering  Dept. 

Texas  A  <fc  M  University 
College  Station,  TX  77S43 

10.  Mr.  Thomas  Momiyama  1 

Naval  Air  Systems  Command 

Aircraft  Division  -  Research  and  Technology 
Air  °31 

Washington.  DC  20360 


11.  Mr.  Marry  Berman 

Naval  Air  Systems  Command 

Aircraft  Division  -  Research  and  Technoloev 

Air  931 

Washington,  DC  20360 

12.  Mr.  Donald  K.  Johnson 
Air  Force  Flight  Test  Center 
6520  Test  Group  ENF 
Edwards,  AFB,  CA  93523 

13.  Prof.  R.  M.  Howard 

Department  of  Aeronautics  and  Astronautics,  Code  67HO 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

14.  LCDR  Robert  W.  Renoud 

USCG  Aircraft  Repair  and  Supply  Center 
Elizabeth  City,  NC  27909-5001 


112 


